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ABSTRACT 



This paper reports an investigation of the dynamic behavior 
of a pin-ended Inelastic column subjected to a half-sine dis- 
placement pulse at one end applied with small eccentricity. 

The rates of end displacement are such that the effects of 
axial inertia are considered negligible. The column is re- 
placed by a lumped parameter mathematical model and the equa- 
tions for the model are solved by a high speed digital com- 
puter. The stress-strain relation during loading is assumed 
to be of the Ramberg-Osgood type. Unloading occurs along a 
line with slope equal to that of the tangent at the origin. 
Slenderness ratios of 50, 70, and 100 and two different values 
of eccentricity are considered in the investigation. The max- 
imum load supported by the column and the residual lateral 
deflection at the center of the column are found as functions 
of the displacement pulse height and of the displacement pulse 
duration. 

It is shown that an inelastic column, like the elastic 
column, can support dynamic loads considerably in excess of 
the static load capacity. It is also shown that within the 
range of pulse durations studied, the smallest residual de- 
flections generally accompany the largest maximum loads. 
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INTRODUCTION 



In the study of structures the column ranks as one of 
the more provoking; elements as judged from the amount of 
literature published on its behavior. The development of 
elastic column theory is comparatively advanced and dates 
back to 1757 when 2uler published his analysis of the strength 
of columns. In more recent years Inelastic column theory has 
been developed to the point where a satisfactory explanation 
has been obtained for the phenomenon. Shanley /l/* has been 
credited with revising and reinterpreting inelastic column 
theory from previously accepted erroneous conjectures. Util- 
izing studies of a highly idealized column he concluded that 
th9 tangent modulus load is the maximum load at which an in- 
itially straight, centrally loaded column will remain straight; 
loading beyond this load causes bending which produces per- 
manent deformation. 

Tne dynamic behavior of structures has become an increas- 
ingly important engineering consideration in recent years and, 
with the aid of more sophisticated devices for the problem 
solution, many more phenomena can be Included in an analysis 
than could previously be reasonably considered. The behavior 
of dynamically loaded elastic columns has received consider- 
able attention and a number of investigations have resulted, 

Hoff and associates /2, 3/ considered the case of an elas- 
tic column of uniform cross section, initially curved in the 

* Numbers in / / refer to entries in the bioliography 
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shape of a half-sine wave, subjected to constant velocity 
loading such as that encountered in compression tests in 
commercial testing machines. It was shown that the maximum 
loads supported by the column are greater than the Euler load 
when the loading is rapid, the column is slender and the ini- 
tial deflections are small. Sevln /4/ confirmed the previ- 
ous results while including the effects of axial inertia 
which were not considered by Hoff. He Inferred that, for an 
elastic column, axial inertia effects are of negligible im- 
portance with regard to the gross behavior of conventional 
structural columns, regardless of the initial deflected shape, 
end constraint, or type of axial loading. 

Taylor /5/ investigated the case of an eccentrically 
loaded elastic column subjected to a half-sine force pulse 
and he included the effects of axial inertia. Using a limit- 
ing value of extreme fiber strain as the failure criterion 
he determined the column failure load quantitatively. He 
demonstrated that elastic columns will support loads much 
greater than the Euler load when subjected to rapid loading 
and that axial inertia effects become significant with this 
type of loading as the rapidity of loading is increased. 

Gerard and Becker /6/ studied column behavior under con- 
ditions of impact, by utilizing the properties of a propa- 
gated elastic compressive stress wave produced on failure 
of a tension specimen. They concluded that a column can 
momentarily support a dynamic compressive stress of any 
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magnitude and that buckling may be confined to a length less 
than the entire length of the column. 

The inelastic column was investigated by Chawla /7/ 
in essentially the same manner as Hoff (with an Initial 
curvature or effective eccentricity) and he found that for 
relatively slow rates of end displacement the maximum column 
load was less than the tangent modulus load. Brooks and 
Vv’ilder /8/ studied the effects of dynamic loading on an 
Idealized H-section inelastic column with an initial lateral 
deflection using a constant velocity end displacement. From 
their results they concluded that the static maximum load 
may be employed as a conservative estimate of the maximum 
load of a column regardless of the rate of end displacement. 

Investigations of the dynamic behavior of inelastic 
columns are limited in number and the problem of determin- 
ing the behavior of inelastic columns subjected to pulse 
loading appears to have been neglected. Yvhereas constant 
velocity loading leads to the ultimate failure of the column, 
pulse loading may not necessarily result in failure. Pulse 
loading allov.’3 the column to recover from a deflected posi- 
tion with an amount of residual deflection while retaining 
a capacity to sustain subsequent load. The problem con- 
sidered in this investigation is that of determining the 
dynamic behavior of a pin-ended inelastic column subjected 
to a half-sine displacement pulse at one end applied with 
small eccentricity. The problem involves the prediction 
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of the maximum force developed In the column for a specified 
pulse neight and duration of loading. The problem also In- 
cludes the prediction of the residual lateral deflection 
resulting from the pulse loading. The effects of pulse height 
and duration are investigated for two values of eccentricity 
and three column slenderness ratios. 
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2 • DEVELOPMENT OF_ THE EQUATIONS FOR THE MATHEMATICAL 

MO DEI . 

The column considered in the investigation is initially 
straight, pin-ended, and of constant cross-section,, The ef- 
fects of axial inertia and shear strains are neglected, but 
non-linear axial strain components due to bending at the 
center portion of the column length are considered. The 
effect of a non-linear stress-strain curve on loading and 
departure from this curve on unloading is also included. 

The loading on the column is produced by axial displace- 
ment of the pin at one end. The loaded end is free to rotate 
about the pin, which is eccentrically placed. The other end 
is constrained by a fixed pin of equal eccentricity. The 
end displacement is imposed by a moving ram executing a half- 
sine displacement pulse. The ram can exert only compressive 
force on the column (see Fig . 1). 

The problem is developed and solved, not in terms of the 
real column, but in terms of the lumx>ed parameter mathematical 
model shown in Fig. 1. The equations for the model are then 
solved in a high speed digital computer. The force-deflection 
characteristic of the springs is derived from the material 
stress-strain diagram. The assumed behavior is indicated in 
Fig. 2 where the solid curve OABD represents the variation 
of compressive stress (T under montonically increasing com- 
pressive strain 6 . If, upon reaching point B, the strain- 
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Pig. 1* Geometry of the Mathematical Column Model, 




I 



6 



ir.g is reversed, the path will be along BC whose slope is 
that of the tangent to the solid curve at 0. A second re- 
versal will cause the compressive stress to increase again 
following CB. After again reaching the curve at B, the path 
will be along the solid curve toward D 0 Additional strain 
reversals would result in additional "detours" parallel 
to BC. In the study the stress-strain relationship after 
departure from the compressive stress- strain curve is con- 
sidered linear into the tensile stress region. The maximum 
tensile stresses in the column model are limited to values 
near the nominal yield stress of the material. It is pre- 
sumed that any tensile stress in the column model springs 
is insufficient to cause plastic straining. 

This behavior is incorporated in the springs by taking 
the stiffness of each as K=A, E T /JL , where A, is half the col- 
umn cross-sectional area, E T is the tangent modulus, and 
is the unstrained length of a spring. Again referring to 
Fig. 2, at any point on the solid curve E r is the slope of 
the curve. For points below the curve E~=E (Young's Modulus). 

This model having only two degrees of freedom was cnosen 
for the investigation in order to limit the computer solution 
time to reasonable values, but yet to retain the effects of 
Inelastic column behavior. It is felt that this model is a 
good compromise in that it approaches the deflected shape 
of the inelastic column and also retains the effects of 
lateral inertia and of the longitudinal strains on each side 
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of the column due to both the axial load and bending. 

From the free body diagram of half the column model 
shown in Fig. 3, using d'Alembert’s principle, the follow- 
ing equations may be written: 



P = f; + Fz 



(1) 



Q—1 1 



dt* 



( 2 ) 



P (e cos-©- + sine-)- Q(±*. cos©- - e sin-e) = 

(Fi cos ^ 

Dividing the latter by cose and substituting for Q from 
equation (2) 



P(e 



ytan-0-) - 



dt* 



— e "tan-e-) - 



(3) 
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From kinematic considerations of the column model in 
a deflected position the following equations may be written 



c£, = + 2.^ sin-e 



(*) 



d z = c$m - z sin© 



( 5 ) 



where c£, and are the axial shortening of Spring 1 and 
Spring 2 respectively and c( m is the average axial short- 
ening at the center of the model. Further 



= + 2. • 4f-(j -cos©) + z e sin© 



(6) 



y =-jr sm ' e ~ e(i-cose) 



( 7 ) 



where x is the end displacement at the loaded end of the 
model and y is the lateral deflection at midlength in the 
model. 

The following approximations are consequences of assum 
ing that is small and e/L is small: 

~ » e t a n © 

sin-© » cos©) 

t an-© ^ sin© ^ © 
l - cos © » 4r ~& Z 

An evaluation of the errors resulting from these approxima- 
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tions is given in App. VI. 

Equation (7) now becomes: 

7 = ( 8 ) 

Now, substituting the approximations, together with equa- 
tion (8), into equations (3) through (6) and rearranging, 
results in 

^ = 4E Ce+ y)-(F;-F,)^ (9) 

c£i ■ dm + y 

d - dm - y 

Solving for and substituting in the above expressions 
for c£, and gives 

<£, = x- - -try 1 - y + ^ y 
-ty 1 - ^-y -*s§y 

Since the model has but two degrees of freedom, only a 
limited correspondence may be established between column and 
model parameters. Energy considerations are used for this 
purpose. 



( 10 ) 

( 11 ) 
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' 















To determine this correspondence both the column and the 
model are considered for the elastic case. Figure 4 indicates 



an elastic column in a deflected position during free vibra- 
tion. It has an assumed half-sine deflected shape and the 
motion is considered sinusoidal with time. 

The maximum potential energy is found from the elastic 
strain energy due to bending and is 

,2. 



El 



Potential Energy = 






dx 



where 412 is the curvature of the axis of the column and I 
d x 2 - 

is the moment of inertia of the column cross-sectional area. 
I is equal to ^"A, where A is the column cross-sectional 
area. This can be shown to be 



Potential 



Energy = 31 



4 El 



4 Li 



3 J 



2 . 



(column) 



The maximum kinetic energy of the column, neglecting 
axial inertia effects, is 

Kinetic Energy =-| C 

'o 

which can be shown to be 



Kinetic Energy 



=i 



A (On y 




(column) 



where ^ is the mass per unit length, <jj n is the circular 
frequency of first mode lateral vibration of the column and 
Lis the column length. Setting the maximum potential energy 
equal to the maximum kinetic energy and solving for the circular 
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> ^ 



y 




y 



y m Bln ~ 



Fig. 4. Elastic Column In Free Vibration. 




Fig. 5* Elastic Column Model In Free Vibration. 
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frequency of first mode lateral vibration of the column re- 
sults in 







( column) 



Figure 5 shows the model with non-eccentric pin connec- 
tions in a deflected position. The stiffness of each spring 
1 s and the motion of the mass is again considered sinusoi- 

dal with time. The maximum potential energy as found from 
the elastic strain energy in the springs is 



Potential Energy = % + <£ z ) 



Considering the model to be in free vibration without axial 
load the average spring deflection is zero or c£, = - <^ z • 
Under these conditions it can be shown from equations (10) 
and (11) that 




and the model maximum potential energy becomes 
Potential Energy = 1 . (N ^ — y* 



/m 



(model) 



r \i 

The maximum kinetic energy of the model neglecting axial 
inertia effects is 



Kinetic Energy = 



(model) 



Now, setting the maximum potential energy equal to the max- 
imum kinetic energy and solving for the circular frequency 
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of first mode lateral vibration of the model results in 



To determine the undeflected length of the model 
springs the potential energies of the column and model are 
equated. Taking the center deflections of column and 
model to be equal, Si. is found to be 



To determine the relation between column and model 
masses the kinetic energies of the column and model are 
equated. Again taking the center deflections to be equal, 
m is found to be 




With the two previous correspondences it can be shown that 
the circular frequencies of the column and the model are the 
same. 

Now, considering the differential equation of motion 
(as previously derived for the model) for the elastic case, 
together with the results of the elastic energy development, 
Euler load and end displacement correspondences between the 
column and model are determined. These will be used later in 
the analysis to define non-dimensi oral load and displacement 
parameters. If the acceleration and the eccentricity of the 




(model ) 



* = ^L 



( 12 ) 
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column model are both zero, equation (9) becomes 



Py = (F,-F 4 )^ 



(») 



Equation (14) applies to static buckling of an elastic column 
with zero eccentricity. Now, for the elastic case 

= and 

From equations (10) and (11) 

Thus 




Substituting this into equation (14), the Euler load P e of 
the column model is found to be 



P = JT 4 A, C 
e 4 (L/^) 2 - 

where L/^ is the slenderness ratio of the column. In terms 
of the total cross-sectional area the equation becomes 



p _ JI 4 A Z 

c ‘ s' TL/O 1 



(15) 



The Euler load, in terms of the Euler displacement, is 

P z ~ ZK X e 



and, substituting for K and P E 




the Euler displacement becomes 

(16) 
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With these results determined from elastic behavior of 
the column and model, the equations for the strains and for 
dynamic equilibrium of the inelastic column model are develop- 
ed in a convenient non-dimensional parametric form. From 
equation (10) the strain in Spring 1 is 






= _cL- — nl cL 

l 32. L 



or 



£ =-H , 

3 2. L L 






L3 J 1} J LT 

Defining a non-dimensional lateral deflection as V= y/e and 
making the substitution y=Ve, together with the substitution 
for the Euler displacement from equation (16), the strain may 
be expressed as 
■4 



r — — — — — JL _ _L ( V 

^ I 8 |_ 2 . v ; v V ^ x + ^ J 



(17) 



Applying this same reasoning for the strain in Spring 2 re- 
sults in 



TT 



6Xl/(oV 






J 



(18) 



From these equations for the strains at the center por- 
tion of the column model it can be recognized that the first 
terms are the dimensionless end displacements, the second 
terms are the corrections to this end displacement due to 
the angular rotation of the rigid rods, the third terms are 
the corrections due to the change of the effective eccentri- 
city with lateral deflection, and the fourth terms result 
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from bending 



The relation between the strains and stresses during 
initial loading of the column model is given by the Ramberg- 
Osgood /9/ equation: 



where 6 is the compressive strain. O’ is the compressive 
stress, CT\j is the yield stress of the material and n is 
the shape parameter of the stress-strain curve. By defining 



This is the equation which is used to determine the stress 
in each spring as the model is initially loaded. For the com- 
pressive stress-strain relation considered in this investi- 
gation the shape parameter is n^lO and the dimensionless 
yield stress is Sip0.004. The stress-strain curve is shown 
in Fig. 6 and is typical of 24S-T aluminum alloy material. 

As each of the springs unloads, it leaves the stress- 
strain curve following a straight line of slope 1 . 

As subsequent loading takes place in either spring the same 
straight line is retraced and, as the stress reaches a level 
that puts it on the stress-strain curve again, it follows 
the curve until further unloading occurs. 

The differential equation of motion is developed in 




a non-dimensional stress S as the previous equation may 

c. 

be written as 



€= S 




( 19 ) 
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nondimenslonal form by dividing equation (9) by the eccentri- 
city e and the square of the circular frequency of first mode 
lateral vibration CO n and substituting from equation (1) for 
P. After- rearranging terms the equation becomes 



cT C//e) 
^ d t*- 



TT 






A dimensionless time T is defined as T = (_O n t, so that 

(jJ z (dt) z = (dT) z . Then — U- becomes ^ . 

d t l d 

By making the substitution V for y/e this becomes which 

d T 

♦ * * * 

is symbolized as V. By substituting V and V and recognizing 

p 

the form as the previously defined non-dimensional stress, 

the differential equation of motion becomes 



v = ^ (L/tf [( s, + so(y+ o - (s, - s t ) 



( 20 ) 



The two terms in the brackets can be recognized, respective- 
ly, as proportional to the external bending moment due to the 
force on the end of the column, and the internal bending mom- 
ent due to the spring forces. 

In order to solve this equation it is necessary to specify 
the forcing function. This is accomplished by specifying the 
displacement of the ram (as shown in Fig. 1 ) as a function of 
time. When the ram is in contact with the column model the 
ram displacement is the end displacement of the column. If 
the ram is not in contact with the end of the column, the 
force is zero (S|=-S^) and the column exhibits free vibration. 
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In either case the conditions of compatibility and dynamic 
equilibrium are maintained. 

The ram displacement is a half-eine pulse of height X 
and duration . To enable a selection of dimensionless 
parameters having a relationship to the column the follow- 



ing dimensionless quantities are defined: 



(a) 



the height of the half-sine ram displacement pulse 
as a fraction of the Euler end displacement; 




(b) the ratio of the natural period of first mode 
lateral vibration of the column to the duration 
of the half-sine pulse; 




With these dimensionless quantities and the previously defined 



dimensionless time. The dimensionless ram displacement 
becomes 



~X ~ s \ n <y>T 



and 



for 



O — T 



^ 2JX 
“ $ 



( 21 ) 



"X. = O for T > 

V 

To determine the maximum force and residual lateral 
deflection resulting from a specified ram displacement it is 
necessary to compute the dynamic response of the model to the 
loading. This is accomplished through a numerical integration 
process applied to the differential equation of motion and a 
numerical iteration scheme to satisfy the relations among the 
strains, stresses, end displacement and lateral deflection. 
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3. THE NUMERICAL SOLUTION 



The equations for the model are solved utilizing FOR- 
TRAN programming and a Control Data Corporation 1604 high 
speed digital computer. The investigation is accomplished 
with two basic programs, COL 1 and COL 2. The former is em- 
ployed primarily to obtain graphical output directly to study 
the dynamic response of the model to a specified loading. Pro- 
gram COL 2 is employed to determine the numerical values of 
maximum force and residual lateral deflection for a specified 
loading on the model. The complete programs in FORTRAN lan- 
guage, together with the program notation, are given in Appen- 
dix IV. 

The numerical integration procedure is the same for 
both programs, but the program, COL 2, employed to determine 
maximum values is more comprehensive. For this reason the de- 
tailed discussion is limited to this program. A discussion 
of the differences between these programs follows the detail- 
ed discussion. Figure 7 is a simplified block diagram of 
Program COL 2 and Fig. 8 is a simplified block diagram of the 
subroutine used in both programs to determine the stress from 
the strain. 

Referring to Fig. 7, the input to the program includes 
the following dimensionless parameters: 

(a) L/^ (slenderness ratio) and e/o (eccentricity ratio) 
to prescribe the column model. 

(b) cx. (ram pulse height) and (ratio of natural period 
of first mode lateral vibration to ram pulse duration) 
to prescribe the forcing function on the model. 
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Fig. 7* Simplified Block Diagram of Computer Program COL 2 
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(c) n ( Ramberg-Osgood shape parameter) and Sy (yield 
stress of the material) to prescribe the compressive 
3tress-strain curve, 

(d) AT (time Increment ) to prescribe the Incremental 
values In the numerical Integration process. 

The Initial conditions are then set to start the numerical 

Integration process. 

Henceforth, the dimensional quantities and the dimension- 
less quantities are referred to synonomously except where con- 
fusion may result, i.e., lateral deflection for dimensionless 
lateral deflection, etc. 

The numerical Integration scheme used in the programs 
1 8 a Runge-Kutta process /10/. This process predicts the 
lateral deflection V and the lateral velocity V at a time 
T+AT from four estimates of the lateral acceleration V cal- 
culated in the time increment AT, Accuracy Is of the order 
of ( A T) 3 . Within this predictor loop further calculations 
are made to satisfy the relations among the strains, stresses, 
end displacement and lateral deflection. 

The ram displacement Is computed from equation (21). 

The values of ram displacement and of lateral deflection are 
then used to determine the strain in Spring 1 from equation 
(17). The main program then calls Subroutine STRESS, Because 
this routine is written with dummy variable Inputs to accommo- 
date the strains and stresses In both Spring 1 and Spring 2, 
the following discussion is presented based on a general 
strain € and a general dimensionless stress S. 



% 
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Subroutine STRESS, shown in Fig. 8, is entered with the 
present value of strain £ and the previous values of stress S 
and strain £ a calculated for the preceding time T-AT. The 
solid curve in Fig. 9 represents the dimensionless compres- 
sive stress-strain curve. The dashed vertical lines indicat- 
ed by A, B, and C in the figure are representative of con- 
ditions which must be considered in arriving at a value for 
the stress S from a value of the strain € . 

Referring now to Figs. 8 and 9, if the present value of 
strain is negative this condition is represented by A. The 
spring has unloaded from the stress- strain curve following 
a straight line with a slope equal to unity and the stress 
is calculated from 

S = $ a + (,£ “ € a ) (22) 

This equation is also used to determine the stress S if €.= 0 . 

If £ is positive there are three possibilities to be consider 
ed„ These are: 

(a) loading is occurring along the compressive stress- 
strain curve; 

(b) unloading is occurring on a unity slope line; 

(c) loading is occurring on a unity slope line. 

To determine which of these possibilities exists the incre- 
mental change in strain, A£= 6 -£ a , Is calculated. If 

this quantity is negative this condition is represented by 
B. The spring is unloading on a unity slope line and the 
stress is calculated from equation ( 22 ). If the incremental 
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Fig, 8. Simplified Block Diagram of Computer Program for 
Subroutine STRESS* 




Plfi* 9* Stress-Strain Considerations In Computer Program 
for Subroutine STRESS* 




/ 
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change in ©train is positive two possibilities exist* Either 
the stress is on the compressive stress-strain curve or the 
stress is increasing on a unity slope line* This situation 
is represented by C in Fig* 9* The stress level is deter- 
mined by calculating a value of stress from equation (22) 
and by calculating a value of stress from the Ramberg-Osgood 
equation (19) and then selecting the smaller of these two 
values as the stress existing in the spring* If the incre- 
mental change in strain is zero, the stress S is equal to S 3 , 

To calculate the value of stress from the Ramberg-Osgood 
equation it is necessary to compute the real, positive root 
of the equation* This is accomplished by a Newton-Raphson 
Iteration /10/ where the strain f, is used as the initial 
estimate for the root Sj. * The stress S could have been de- 
termined from the strain € each time by calculating a value 
from equation (22) and a value from the Ramberg-Osgood equa- 
tion (19) and selecting the smaller value* However, this 
would involve an iteration for each stress calculation, thus 
increasing the computer time* Therefore, this method was not 
used* 

After S is determined the value for the strain £ z is 
calculated from equation (18)* The main program again calls 
Subroutine STRESS to determine the value of S* „ 

At this point the program calculates the sum S,+S z , 
which is directly proportional to the force on the end of 
the model, to determine if it is compressive (positive) or 
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tensile (negative). As previously discussed for the model 
configuration (Fig. 1), the ram cannot transmit a tensile 
force to the model. Therefore, if this sum is negative the 
program Initiates an iterative Interpolation scheme to deter- 
mine the end displacement of the model which will cause the 
sum Si+Sj, to be zero. This involves recalculating the strains 
and then the stresses for each end displacement until the force 
is zero. The model is then in a state of free vibration. This 
free vibration continues until the ram again contacts the column 
or until the calculation terminates. 

Once these states of stress (end force either positive 
or zero) are determined, the lateral acceleration V is cal- 
culated from equation (20). After four passes are made to 
this point in the program to determine the values for the 
Runge-Kutta predictions, as previously discussed, the lateral 
velocity V and the lateral deflection V are calculated 
for the end of the time interval A T. 

At this point the program computes the force which can 
be shown to be 



and a continuous check is made to determine the maximum value 
of this quantity which is stored for print-out. This contin- 
uous check is also made to determine the maximum value of 
compressive stress in each of the springs. These values are 
stored for use in strain energy calculations to be made later. 




2 ? 



A continuous check Is also made to determine the extreme stress 
(tensile or compressive) In each spring for prlnt-out. The pro 
gram then calculates the ram work done on the model for the pre 
ceding time Increment and adds this quantity to the previous 
work done by the ram. 

This whole process continues until the ram displacement 
pulse has terminated and the column Is left In a state of 
free vibration. Then a continuous check Is Initiated to 
determine the time at which S|=S2.=0. The lateral deflec- 
tion at this point Is the residual deflection y>- In the 
model, as the lateral acceleration, the longitudinal end 
force and the Internal moment are all zero, and the residual 
strains determine the permanent deformation that has resulted. 

After the residual deflection has been calculated a 
switch is set to detour this portion of the program. The 
new path Includes a continuous check to determine the time 
at which the lateral velocity is zero (point of extreme lateral 
deflection and zero kinetic energy). At this point the net 
energy Input to each of the springs is calculated from the 
area under the stress-strain curve using the present values 
of stress and the previously stored maximum values of the com- 
pressive stress. This energy is a combination of the strain 
energy of plastic deformation and of the potential energy due 
to existing elastic strain in the model. This energy should 
be equal to the ram work done on the model. 

When the energy has been calculated the program then 
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prints all input quantities and the results. The results 
include the values of P^/Pj. ,y*-/^ ,S lrnax , S ima% and the 
time T that each of these occur, as v.'ell as the work and 
energy. The times are included to complement the time plots 
from the other program as an aid in visualizing some of the 
dynamic behavior and to determine the problem time termina- 
tion of the program. 

The program employed to obtain graphical output directly 
is Program COL 1 (the FORTRAN language listing is given in 
Appendix IV) . This program differs basically from Program 
COL 2 described in the preceding paragraphs in that it stores 
the variable quantities for each time T in array form instead 
of storing only the maximum values. These array quantities 
are then available for Subroutine GRAPH (not shown in the 
listing but available on the FORTRAN Library Tape at the U.S. 
Naval Postgraduate School, and also available from the Control 
Data Corporation). Subroutine GRAPH is designed to accept 
a maximum of 900 points for each variable to be plotted. Thus 
the program is written to calculate only 900 points while main- 
taining the necessary accuracy for the numerical integration 
process compatible with the time increment AT. 

In addition to the input data required in Program COL 2, 
axis scaling and multiple curve data, graph title data and 
curve label data are required for Frogram COL 1. In addi- 
tion to the graphical output the program prints the input 
quantities and the values of T, x/x £ , , S, , S z , y/^ 

and F/p E for every fifteenth point. As all of these quantities 
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have been stored as 9G0 point arrays any combination of these 
variables may be plotted as a function of another. Subroutine 
GRAPH converts these 900 point arrays to input data on a mag- 
netic tape for a Control Data Corporation 16G computer which 
is programmed to provide output as input to a plotter. 

An example of one set of curves used in the investi- 
gation is shown in Appendix I. The example includes three 
graphs with two curves on each graph. The first graph is a 
plot of €, and versus T, The second graph is a plot of 
S t and versus T and the third graph is a plot of y/^ and 
F/P e versus T, Other examples of graphical output are shown 
in Appendices II and III. However, the graphs shown in 
Appendix II, which are for very slow loading approaching the 
static case, have been obtained with a modified version of the 
program. The modification Includes calculating 1800 points 
to maintain the accuracy for the numerical integration pro- 
cess, then re-storing every other point of the 1800 point 
arrays as 900 point arrays to accomodate Subroutine GRAPH. 
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RESULTS AND DISCUSSION 



a. Character of the response 

An understanding of the results is dependent upon an 
appreciation of the distinctive features of the loading and 
their effects upon the response. Loading Is accomplished by 
a forced displacement of one end of the column in the form 
of a half-sine pulse. This pulse is characterized by the 
dimensionless parameters cx. and „ The quantity is the 
ratio of the displacement pulse height to that end displace- 
ment which, in the absence of lateral deflection, corresponds 
to the Euler buckling load. The parameter is the ratio of 
the first mode natural period of lateral vibration of the col 
umn model to the pulse duration. Thus, ex measures the ampli 
tude of the loading and measures its rapidity. 

Consider the very slow loading (Q}«l) of a slender 
column to a maximum end displacement beyond the Euler crlti- 
al value (oo I ). At axial loads above the Euler limit 
the straight form of the column is unstable. To initiate 
lateral deflection of the column model the end pins are 
eccentrically located. The effect of this eccentricity is 
to produce bending at loads below the stability limit of the 
column. The maximum force (P rn ) developed in the column 
becomes less than the Euler load (P t ). However, for very 
slow loading the maximum force developed in the column model 
is the static maximum load. Increasing values of CX above 
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that corresponding to the Euler critical value do not change 
the value of the maximum force P^ . 

For increasing loading rates (increasing values of ) 
the lateral motion of the column is retarded by the inertia 
of the mass and any meaningful correspondence to static be- 
havior is lost. However, the maximum force developed in the 
column is related to the rate of loading. Also, under pulse 
loading the inelastic column recovers from a deflected posi- 
tion with an amount of residual deflection and retains a 
capacity to sustain further load. This capacity may be 
measured in terms of the amount of residual deflection 
which results from the loading. This residual deflection is 
equivalent to additional eccentricity for any subsequent load- 
ing. The maximum force P^ and the residual deflection y r appear 
to be the most important performance criteria for the dynamic 
behavior of the inelastic column subjected to a pulse loading. 
Dimensionless forms of these two quantities are utilized in 
the graphical form of the results. 

b. Range of variables 

Solutions for Pm/Pe. (maximum force) and for yr/^ (residual 
deflection) have been obtained for a range of from 0.01 to 
1.0 and for a range of cx. dependent on the column slenderness 
ratio. The results include solutions for three values of col- 
umn slenderness ratio (L/^ ) and two values of eccentricity 
ratio (e/^ ). For slenderness ratios of 70 and 50 the mini- 
mum pulse height considered was c<= 1. For a slenderness 
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ratio of 100 the minimum pulse height considered was cx= 2. 

For values of cx less than these, inelastic action is insig- 
nificant, The maximum values of oc. considered are based on 
the maximum value of tensile stress occurring in the column 
model springs. For the investigation the stress-strain re- 
lationship after departure from the compressive stress-strain 
curve is considered linear into the tensile stress region. 

For this reason the maximum tensile stresses in the column 
are limited to values near the nominal yield stress of the 
column material. For slenderness ratios of 70 and 50 the max- 
imum pulse height considered was cx. = 2,5, For a slenderness 
ratio of 100 the maximum pulse height considered was cx=* 4. 

For values of cx greater than these the maximum value of ten- 
sile stress occurring in the column model is considerably in 
excess of the nominal yield stress for some rates of loading. 

c. Maximum force 

Figure 10 is a plot of families of curves of maximum force 
P lri /P E versus P> for a constant pulse height cx . The families 



Include slenderness ratios of 100, 70, and 50 for an eccentri- 
city ratio of 0.05. Figure 11 is a plot of the same families 
of curves for an eccentricity ratio of 0,01, The lover limit 



case. The upper limit of Q ~ 1.0 is approaching the region 
where axial inertia effects become important in the elastic 
column as reported by Taylor /5/« As the inelastic column 




is for slow loading rates approaching the static 



model does not include the effects of axial inertia 
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1 8 the upper limit considered in this study. 

Figures 10 and 11 show that as the rate of loading is 
increased the maximum force developed in the column increases. 
This is due to the delay in the development of lateral deflec- 
tion because of lateral inertia. These figures also indicate 
that as the rate of loading is decreased the maximum force de- 
veloped in the column approaches a limit asymptotically. In 
Fig, 10 this limit may be observed for each slenderness ratio 
as all of the curves converge to a single value at Q = 0.01. 
This limit is the static load capacity of the column for an 
eccentricity ratio of 0.05. In Fig. 11 this limit is not 
attained. The spread of points indicates that for this eccen- 
tricity (e/(> = 0.01) there is considerable dynamic effect at 
Q = 0.01. However, the curves of minimum cx. f 0 r each slen- 
derness ratio at = 0,01 appear to approach a limit of the 
static load capacity of the column fcr an eccentricity ratio 
of 0.01. A comparison of these figures confirms that as the 
eccentricity increases the maximum force developed in the 
column decreases. 

d. Slow loading 

It is possible to make some observations concerning 
the values obtained at the low end of the range by con- 
sidering the static case of a centrally loaded inelastic 
column. For this ca3e the maximum load the column can sup- 
port is between the tangent modulus load P T and the reduced 
modulus load P^ . In terms of the Euler load the dimensionless 
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tangent modulus load Is 




( 23 ) 



where E T is the tangent modulus or the local slope of the 
compressive stress-strain diagram. The tangent modulus de- 
fines the value of stress at which bending begins for the 
static case of a centrally loaded Inelastic column. In 
terms of the model constants developed previously this 
critical stress is 



0”cr — 



_ TT 4 E- 



8 ( L _/^) 4 

or in terms of dimensionless stress this becomes 



Sr - -Hi ( 1— ' \ 

0cr ~ ft T \L/^i 



(24) 



The ratio of Young's modulus to the tangent modulus in terms 
of the Bamberg- Osgood parameters is 






( 25 ) 



where n = 10 (the shape parameter) and = C.0C4 (the dimen- 
sionless compressive yield stress) are the parameters used 
in the investigation. By choosing a value of S cy - and sub- 
stituting this for S in equation (25) the ratio E/E T is de- 
termined. By substituting the values of this ratio and S C r 
in equation (24) the slenderness ratio is determined. In 
this way a value of the ratio Ej/E is determined correspond- 
ing to each of the slenderness ratios 100, 70, and 50. This 
value of Sr/S equals P T /P e from equation ( 23 ). 
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In terms of the Euler load the dimensionless reduced 



modulus load is 



P e £ 



(26) 



where E R is the reduced modulus. The reduced modulus de- 
fines a maximum value of stress which can be approached but 
not exceeded for the static case of a centrally loaded inel- 
astic column. In terms of the model constants developed pre- 
viously this critical stress is 

rr - t 4 

^ " 8 

or in terms of dimensionless stress this becomes 



S, r . = Oil Sl/jn 8 

° cr 8 E lL /^/ 



(27) 



The reduced modulus 13 dependent on the tangent modulus and 
the geometry of the column. For the column model it can 
be shown that the reduced modulus is 
r - 

t ' K " E + E t 

or In dimensionless form this becomes 



Ea. z Er/F. 

E \ + E t /E 



(28) 



and thl3 value is equal to P R /P e from equation (26). By 
choosing a value of S c> - and substituting this for 3 in equa- 
tion (25) the ratio E t /E Is determined. The ratio E t /E also 
determines E*/E from equation (28). By substituting the values 
of E r /S and 3 Cr in equation (27) the slenderness ratio is deter- 
mined. In this way a value of the ratio E R /S is determined 
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corresponding to each of the slenderness ratios 100, 70, and 
50. This value of E^/E equals P*/P E from equation (26). 

The values of the tangent modulus load and the reduced 
modulus load determined from these two methods are: 

for L/^ = 100, P T /P £ = 1.000 and P R /P e = 1.000; 

for L/^ = 70, P T /P t = 0.969 and P R /P E = 0.977; 

for L/<^ = 50, Pr/Pf- = 0.649 and P R /P p _ = O. 69 I . 

From Figs. 10 and 11 It is seen that In most cases the max- 
imum values at Q = 0.01 are less than the tangent modulus 
load although the spread of points for an eccentricity ratio 
of 0.01 indicates that dynamic effects are more pronounced 
the smaller the eccentricity. For tills reason some of tne 
maximum values are greater than the reduced modulus load. 

To observe further the effect of eccentricity on the 
load capacity of the column at very slow loading rates, a 
study was made to determine tbs behavior of the column under 
a ramp displacement of the end at a loading rate correspond- 
ing to that available to commercial testing machines. In 
terms of the parameters previously used to describe the pulse 
loading of the column the ramp displacement of the head is 

X = Q> T 

where cx. = 1.0 and = 0.01. Computer solution results in- 
cluding a graphical presentation of Strain 1 ( £, ) and Strain 
2 ( 6^ ) versus force (P/P e ) were obtained for the three column 
slenderness ratios of 100, 70, and 50 for eccentricity ratios 
of C.05, 0.01, and 0.001. The graphical results shown in App.II 
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exhibit a general resemblance to the experimental strain 
versus force curves of actual column test data displayed by 
Shanley /l/. From the results given in Table 1 it is seen 
that for eccentricity ratios of C.05 and 0.01 the maximum 
load is less than the tangent modulus load for all three 
slenderness ratios considered. However, the maximum load 
increases with a decrease in the eccentricity ratio. As 
the eccentricity ratio is decreased to 0.001 the maximum 
load is greater than the tangent modulus load but less than 
the reduced modulus load for a slenderness ratio of 50. For 
a slenderness ratio of 70 there is considerable dynamic effect 
and the maximum force exceeds the reduced modulus load by a 



Table I Computer and Theoretical Static Load Comparison 
(For Ramp Loading of Column Model) 





L/<^ = 100 ' 


L/^> = 70 


L/^> = 50 


Static 

Pr/P E 


1.000 


0.969 


0.649 


Static 


1.000 


0.977 


0.691 


Dynamic 
( Computer) 


for 

e/^ = 0.05 


0.952 


0.839 


0.584 


for 

e/<^ = 0.01 


0.983 


0.916 


0.633 


for 

e/(> = 0.001 


1.000 


0.981 


c .665 
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small amount. For a slenderness ratio of 100 there is consider- 
able dynamic effect also but the solution oscillates about an 
average maximum value of P^/ P E = 1,0. 

The tabulated values of P m /P E for an eccentricity ratio 
of 0.05 for a ramp displacement loading are in agreement with 
the maximum force values P m /P E for pulse loading at Q> = 0.01 
in Fig. 10. The strain ( £, and £ 4 ) versus force (P/P E ) curves 
for the ramp displacement loading for an eccentricity ratio 
of 0.05 shown in App. II display no dynamic behavior. In Fig. 

11 (e/(D = 0.01), where the dynamic effects are more pronounced 

at Q> * 0.01 for pulse loading,the forces on the C*.= l curves 
for slenderness ratios of 70 and 50 agree with these tabulated 
results. For a slenderness ratio of 100 the smallest value of 
c<» considered for pulse loading was 2 and this curve appears 
to approach a value P^/P^ = 0.983 asymptotically. The strain 
versus force curves for the ramp displacement loading for an 
eccentricity ratio of 0.01 shown In App. II display an insignif- 
icant amount of dynamic behavior which does not affect the 
values of maximum force P_/P.. . Therefore, these values of 

P^./Pj. may be considered the static load capacity of the inel- 
astic column for eccentricity ratios of 0.05 and 0.01. 

Appendix III contains a compilation of the dynamic re- 
sponse curves stress (Si and 3? ) versus time (T) for a range 
of from C.l to 1.0 and for an eccentricity ratio e/c^ = 0.01. 
This graphical compilation includes the solutions for the min- 
imum and maximum values of pulse height consldex^ed in the 
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study for all three slenderness ratios 0 Included with the 
compilation Is a discussion of the more prominent features 
of the curves,. 



e. Residual deflection 

Figures 12 , 13, and 14 are famlles of curves showing 



height ex. and for e/^ = 0 o 01 and 0 o 05o Figure 12 Is the 
family for L/C = 100, Fig c 13 Is the family for L/^ = 70, 



the values of residual deflection Increase with the rapidity 
of loading to a maximum value then decrease abruptly „ This 
abrupt decrease is because of the delay in the development of 
lateral deflection due to lateral Inert la Q From Fig. 13 
for L/x^ = 70 and Fig. 14 for L/c^ = 50 it is shown that in gen- 
eral the values of residual deflection are essentially con- 
stant as the rapidity of loading increases until an abrupt 
decrease occurs due to the effect of lateral inert la 0 The 
exception for cx = l and L/£ = 70 is due to the small amount 

of Inelastic action that takes place for this pulse height 
at all values of 0} „ The amount of residual deflection may 
be expressed in terms of the amount of inelastic action that 
occurs. From the graphical compilation of stress versus time 



residual 




and Flg 0 14 Is the family 




The points shown 



on these curves are from the identical computer solution 
results as the points shown in Figs 0 10 and 11 » 



From 




shown that in general 



curves in App c Til and the data complied in App 0 V the amount 







Fig. 12. Values of residual deflection yy/^ vs. 

(cx. constant, e/^ °»0.01 and 0.05) for a 
column slenderness ratio L/^ of 100 
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Fig. 13. Values of residual deflection y*-/P vs. Q 
(cx.constant, e/^* 0.01 and 0.05) for a 
column slenderness ratio L/^ of 70. 
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$ 



Fig. 14. Values of residual deflection yr/$ vs. ^ 
( ck. constant , e/<£= 0.01 and 0.05; for a 
column slenderness ratio L/^ of 50. 



of Inelastic action may be determined 0 For a slenderness 
ratio of 100, inelastic action in the column is limited to 
maximum compressive stress levels below the value of the 
yield stress for all solutions considered in the investiga- 
tion. For a slenderness ratio of 50, inelastic action, in 
general, occurs at much higher maximum compressive stress 
level s-near or above the yield stress. For a slenderness 
ratio of 70, inelastic action occurs at maximum compressive 
stress levels Intermediate to those for slenderness ratios 
of 100 and 50. 

A. comparison of Fig. 12, 13, and 14 Indicates that as 
the slenderness ratio decreases the effect of eccentricity 
on the values of residual deflection is less pronounced. 
However, decreasing eccentricity causes a shift of the re- 
sidual deflection curves to longer pulse durations (smaller 
values of § ). This is because dynamic behavior is more 

pronounced the smaller the eccentricity and lateral Inertia 
effects become more important at slower loading rates. It 
may also be observed from these residual deflection curves 
that the maximum values of residual deflection do not necessa- 
rily occur at the slowest loading rates. The minimum values 
of residual deflection generally occur at the fastest loading 
rates. The only exceptions are for the columns of largest 



follow the same trend but are limited by the maximum value of 
the parameter^ = 1.0( shortest pulse duration considered in the 
investigation) . 



slenderness ratio 




and largest eccentricity 




However, the residual deflections appear to 
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CONCLUSIONS 



From the results discussed in the preceding section, 
the following conclusions may be drawn for the rates of 
loading considered in the investigation: 

(a) The inelastic column model, together with a high 
speed digital computer and the peripheral plotting 
equipment, provides a sound and reasonable method 
for the study of the dynamic behavior of inelastic 
columns. 

(b) For very slow loading (long pulse duration) the model 
exhibits behavior in agreement with Shanley's postu- 
lates. 

(c) As eccentricity is reduced, the lateral response is 
slowed. Accordingly, columns with very small eccen- 
tricity show significant dynamic effects at quite 
slow loading rates. 

(d) Within the range of pulse durations studied, the 
maximum load Increases monotonically with decreas- 
ing pulse duration. 

(e) The ratio of maximum dynamic load to static load 
capacity decreases as the slenderness ratio is 
reduced. 

(f) Within the range of pulse durations studied, the 
smallest residual deflections generally accompany 
the shortest pulse durations. 
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APPENDIX I 



Example of Graphical Output from Program COL 1 
Graphical Notation 

The computed points on the graphs are the values of the 
dimensionless variables developed In the basic equations. 

The graphical notation Is as follows: 



Symbol on Graph 



Corresponding Symbols In 
the basic equations 



SLENDERNESS 


L A; 


ECC 




e /'< 


ALPHA 




cx. 


BETA 




§ 


STRAIN 


1 


C 


STRAIN 


2 




STRESS 


1 


s 


STRESS 


2 


s a 


DEFLECTION 


y/^ 


FORCE 




p/p E 


TIME 




T 



Each of the curves on the graph Is labeled to correspond 
to the terms In the graph title. The x-axis Is the axis 
of abscissas and the y-axis Is the axis of ordinates. The 
numbers specifying the coordinates are integer numbers locat- 
ed to the right of the coordinate on the x-axls and above the 
coordinate on the y-axis „ The scale quantity Is the value of 
the first marked coordinate from the origin in either the 
x-direction or the y -direct ion expressed in powers of ten 
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to locate the decimal point . The symbol E In the scale value 
indicates that the number following is the exponential power 
of ten to be applied to the basic scale value, i.e., 1.00E+C1 
is 1.00 x 10 +l . 

Discussion of the curves 

On pages 1-5 to 1-7 following this discussion are shown 
curves of 6, » , S, , S z , y/^ and P/P E versus T. These 

curves are the solutions for the inelastic column model dy- 
namic response to a ram pulse loading with the following par- 
ameters : 

(a) slenderness ratio, L/^ = 100 

(b) eccentricity ratio, e/^ = 0.05 

(c) ram pulse height, cx. = 4.0 

(d) ram pulse duration in terms of = C.2, 

T* = 2 TT/Q> = 10T T 

The discussion includes some of the more prominent features 
of the dynamic behavior of the column model under ram pulse 
loading. 

A comparison of the strain versus time curves with the 
stress versus time curves indicates the amount of inelastic 
action that takes place in Spring 1 and Spring 2. In Spring 1 
the compressive stress level increases above that in Spring 2 
as the amount of bending increases. The internal bending 
moment in the column model is directly proportional to 
S, - Sg_ . As the amount of bending increases Spring 2 de- 
parts from the compressive stress-strain curve along a straight 
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line with unit elope and continues into the tensile stress 
region. The values of stress in Spring 2 remain on this 
straight line for the remainder of the solution. The com- 
pressive stress level in Spring 1 continues to increase to a 
point where it also departs from the compressive stress-strain 
curve along a straight line as the internal bending moment in- 
creases and the force decreases to zero. This force is dir- 
ectly proportional to S + S 0 When the force becomes zero 
the end of the column has departed from the ram which exhibits 
a half-sine pulse displacement. The column model is then in 
a state of free vibration until the end again makes contact 
with the ram. The compressive stress level in Spring 1 again 
returns to values on the compressive stress-strain curve. 

This alternate departure and contact between the column and 
the ram continues until the ram displacement pulse nears its 
termination. The column then continues in a state of free 
vibration where S, = - S i0 The amount of residual deflec- 
tion may then be determined from the mean value of the deflec- 
tion or from the deflection at the time when S = S- = 0, At 
this point the internal forces and the Internal bending moment 
in the column model are zero. 

The effect of eccentricity on the dynamic behavior of 
the column model with all other parameters remaining the 
same may be observed by comparing the stress versus time 
curves in this appendix with the stress versus time curves 
on page III-13 where the eccentricity ratio is 0,01, It may 
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be noted that the essential differences In these curves are 
that the maximum compressive stress levels attained In Spring 1 
and Spring 2 are greater for the smaller eccentricity „ Also, 
contact between the ram and the column terminates sooner for 
the smaller eccentricity 0 
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APPENDIX II 



Theoretical _ and Computer Solution Graphical Comparisons of 
Maximum Loads for the Static Case 

The notation for the graphs displayed in this appendix 
is explained In App c I c 

The graphB on pages II-3 to 11=11 folowing this dis- 
cussion are the curves ? and • , versus P/P for the three 
slenderness ratios considered in the study , As previously 
mentioned, the loading for the study in this appendix is a 
ramp displacement of the end of the column model at a dis- 
placement rate comparable to that in commercial testing ma- 
chines,, The study includes eccentricity ratios of 0,05, 
0,01, and 0 o 00l for each of the slenderness ratios,, The 
tabulated results of this study are shown in Table I in 
Section 4, 

In addition to the previous comments made concerning 
the results, it may also be seen from the graphs that for 
all slenderness and eccentricity ratios considered an amount 
of bending takes place before strain reversal occurs in 
Spring 2 and before the maximum load is reached. This was 
essentially Shanley 0 s point in proving the fallacy of the 
reduced-modulus theory for column action. As stated by 
Shanley /!/ 8 

,,,ocothe reduced-modulus theory is not correct 
for predicting the maximum load up to which a per- 
fect column will remain straight. This is because 



II-l 



it is possible for the column to bend simultaneously 
with increasing axial load a Under such conditions 
it is possible to have bending without introducing 
any strain reversal, upon which the reduced-modulus 
theory depends 0 
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APPENDIX III 



Column Model Dyna m ic Response Cu r ve s 

The notation for the graphs displayed in this appendix 
is explained in App. I. 

The graphical compilation on pages III-3 to III-26 
following this discussion gives the cur ves S t and S z versus 
T for Q = 0.1, 0.2, 0.5 and 1.0 and for an eccentricity ratio 
e/^ = 0.01. Included in the compilation are the solutions 

for the minimum and maximum values of ram pulse height cx. con- 
sidered in the investigation for all three slenderness ratios. 
The graphical solutions are grouped by slenderness ratio with 
the solution for the maximum value of cx. followed by the solu- 
tion for the minimum value of cx» for a single value. The 
order of the (^> values is C.l, 0.2, 0.5 and 1.0, that is, from 
the longest pulse duration of loading to the shortest pulse 
duration of loading considered in the compilation. 

A discussion of some of the prominent features of the 
dynamic behavior of the column model under ram pulse loading 
is included in App. I where the curves of 6, , 6^, S, , S z , y/^ 
and P./P E versus T are displayed for one computer solution. 

The following Is a discussion of some of the more prominent 
features of the dynamic behavior of the column model in terms 
of the stress curve compilation in this appendix. 

In general the amount of inelastic action in the column 
model is Inversely related to the slenderness ratio. For a 
slenderness ratio of 100 inelastic action is limited to cora- 
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pressive stress levels below the value of the yield stress. 
For a slenderness ratio of 50 the maximum conpresslve stress 
levels are near or above the yl9ld stress. For a slenderness 
ratio of 70 inelastic action occurs at maximum compressive 
stress levels intermediate to those for slenderness ratios 
of 50 and 100. 

It may also be noted that, in general, contact is lost 
between the ram and the end of the column (when S ( = - S^) 
before the loading pulse nears its termination for values 
of 0.2 and 0.5 and slenderness ratios of 70 and 50. For a 
slenderness ratio of 100 this occurs at values of (^> = 0.1, 

0 .2 and 0.5. 

It may be seen that the value of maximum tensile stress 
in Spring 2 exceeded the value of the nominal yield stress 
(3y = C.CC4) for the values of 0.2 and 0.5 and ex =» 2.5 
for a slenderness ratio of 50. 
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APPE NDIX _IV 

Comx'uter programs and Not ation 
a. General 

The FORTRAN language computer programs, together with 
the notation, and the program operating instructions for 
both Program COL 1 and Program COL 2 are included in the 
following pages. The discussion and a listing of Program 
COL 1 are Included on pages IV - 2 to IV -10. The discuss- 
ion and a listing of Program COL 2 are included on pages 
IV -11 to IV -20. The FORTRAN language format for each of 
the programs is in standard FORTRAN language to be used with 
a standard FORTRAN compiling routine. The program language 
is compatible with any of the large machine FORTRAN compiling 
systems including the I3M 7090 and the CDC 1604 systems. Each 
line in the listing is toe information contained on one 80 col- 
umn IBM card where columns 1-5 contain any statement number, 
column 6 is used to indicate a continuation card, columns 
7-72 contain the FORTRAN statement, and columns 73-80 are 
used for card identification only. Although the notation is 
basically the same, the subscript notation is different for 
each program. Program COL 1 stores a number of variable 
quantities in array form for print-out and to accommodate 
Subroutine GRAPH, whereas. Program CCL 2 stores only one 
value of tne variable quantities for print-out. 
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b e Notation for Program COL 1 

In the notation for Program COL 1 the subscript J re- 
fers to quantities within the Runge-Kutta prediction loop 
which are calculated for intermediate times within the time 



Increment 


AT. The subscript I 


refers to the 


quant it les 


corresponding to each time T(I). 


The notation 


for Program 


COL 1 is 


as follows: 








Computer program notation for Program 


_C0L_1 


Computer 

Program 

Term 


Correspond ing 
terms in table 
of symbols 


Computer 
Program 
term _ 


Corresponding 
terms in tabl- 
of symbol 3 


A 


n 


SI GY 


S, 


ALPHA 


CX. 


SIG1( J) 


S, 


BETA 




SIG2( J) 




DDV ( J) 


V 


SLNESS 


L/^ 


DEFL(I) 


V 


STRNl(I) 


€, 


DELT 


AT 


ST RN2( I) 




DISPL(I) 


x/x E 


ST RSI ( I) 


3 1 


ECC 


e A> 


ST RS2 ( I ) 


Sz 


EPSI1( J) 




T(r) 


r 


EPSI2( J) 


£-Z 


VSL 


V 


FORCE ( I ) 


P/P E 


X(J) 


x/x E 


o. Operating instructions for Program COL 1 




The 


output from Program COL 


1 includes the 


print-out 



and the array quantities necessary for Subroutine GRAPH. 
The listing for Program CCL 1 does not Include Subroutine 
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GRAPH. However, a listing of Subroutine GRAPH and instruc- 
tions for its use may be obtained from the Control Data Cor- 
poration. To use the program to obtain the print-out quanti- 
ties without the graphical output remove the cards CALL GRAPH. 
To operate Program CGL 1 without Subroutine GRAPH use the cards 
including Subroutine STRESS in the order shown in the listing 
followed by one blank card, then the input data card. 

The read-write statements in the listing are written 
for magnetic tape input and output as follows: 

(1) Input (READ INPUT TAPE 2) 

The input to the program is contained on one card (731C.3) 
and includes the following seven quantities in order: 

I nput quantity Format 



SLNESS (L/q> ) 


E10.3 


ALPHA ( c< ) 


E10.3 


ECC (e/q>) 


E10»3 


BETA ( qb ) 


£10.3 


A (n) 


E10.3 


SIGY (3y) 


E1C.3 


DELT (AT) 


E10.3 



If Subroutine GRAPH is used axis scaling and multiple 
curve data, graph title data and curve label data in 
addition to the above are also needed in the input. 

(2) Output (WRITE OUTPUT TAPE 4) 

The output Includes the input quantities and the vari- 
able quantities for every fifteenth point as follow: 
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Output Quantity 


Format 


SLNESS (lA> ) 


F5 o 0 


ALPHA ( ) 


F5.2 


ECC (e/^ ) 


F5.3 


BETA ( ) 


F4 e 2 


A (n) 


F5d 


SIGY (Sy) 


F6 0 3 


DSLT (AT) 


F5..3 


T 


F5.2 


DISPL (x/x t ) 


E14 0 4 


STRN1 ( f s ) 


214 0 4 


STRN2 (£.) 


E14 0 4 


ST RSI (S, ) 


E14 0 4 


3TRS2 (S z ) 


214 A 


DEFL (y/^ ) 


314 c 4 


FORCE (P/Pc) 


El 4 0 4 
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PROGRAM COL 1 



C PROGRAM TC DETERMINE INELASTIC COLUMN RESPONSE TO HALF-SINE 
C RAM DISPLACEMENT PULSE APPLIEC AT END WITH SMALL ECCENTRICITY 
C TO INITIALLY STRAIGHT COLUMN. WHEN RAM CAUSES TENSION COLUMN 
C EXHIBITS FREE VIBRATION. 

OD I MENS ION T ( 900 ) , C I SPL ( 90 C) , C EFL ( 900 ) , STRN 1 ( 900 ) ,STRN2 ( 9C0), 
1 STRS1 ( 90C ) ,STRS2( 900 , X< 4 ) , QCV ( 4 ) , EPS II ( 4 ) , EP S 12 (4 *,SIG1 (4) , 
2 S I G2 ( 4 ) ,AK( 4,2) ,C (L ) f FORCE (900) 

COMMON A , S I G Y 

READ INPUT TAPE 2 , 5 , SLNE SS , ALPH A , ECC , BET A, A , S IGY , D ELT 
S F0RMAT(7E1C.3) 

T ( 1 ) = C . C 
OEFL ( I )=0 .0 
VEL=0 . C 
STRN 1 (0=0.0 
STRN2 (0=0.0 
STRS1 (C ) =0.0 
STRS2 (0=0. C 
C(1 )=C.O 
C ( 2 ) =C. 5 
C ( 3 ) = C . 5 
C ( 4 ) = 1 . C 
KTENS= 1 
KON T ACT= 1 
M = 0 



l'/-5 



1 = 1 

40 CO 70 J = 1 f 4 

TC = T ( I 1+C( J )«DELT 

VC = DE FL ( I )+C( J)*AK( J-l, 1 ) 

IF( 0. 5*BET A «TC-3. 14 1 59 3) SO, 55, 55 
50 X( J)=ALPHA*SINF(0.5*BETA*TC) 

GO TO 60 
55 X(J)=C.O 

600EPSI 1 C J) = 12. 1761/ ( SLNESS*«2 ) * ( X ( J ) -0. 5* C ECC** 2 )* VC **2* 
1 EtC*VC*( 1 .0-ECC)) 

I F ( J- 1)62,62,63 

62 CALL STRESS(EPStl (J), STRN 1(1-1), STRS1 (1-1), SIGI(J) ) 

STRN1 ( I ) =E PS I 1 (1) 

STRS1 ( I ) =S I G 1 (1 ) 

GO TO 64 

63 CALL STRESS(EPSI 1 ( J ) ,STRN1 ( l ) ,STRS1 ( I ) , SIG 1 ( J ) ) 
640EPSI2(J)=12. 1761/ ( SLNESS**2 )*(X(J)-0.5*(ECC*«2)*VC**2- 

1ECC*VC*( 1 • 0+ECC ) ) 

IF(J- 1)65,65,66 

65 CALL STRESS(EPSI2(J),STRN2( I-1) »STRS2( I— 1 ) , S IG2 ( J) ) 

STRN2 ( I ) = EPSI 2(1) 

STRS2 ( I ) = S I G2 ( 1 ) 

GO TO 661 

66 CALL STRESS( EPSI2 ( J) ,STRN2( I ) ,STRS2( I ),SIG2( J ) ) 

661 IFtKTENS- 1 )67 ,665 ,665 
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665 IFIK0NTACT-1 167 »6 94 » 694 



67 I F ( M- 1 1675,68,685 

i 

675 XI =X( J ) 

P1 = SIC-1 ( J ) +S IG2 ( J ) 

X( J)=CISPL{ 1-1 ) 

M=M+1 
GO TC 60 

68 X2*X(J) 

P2*S I G 1 ( J 1 +S I G2 1 J ) 

X( JJ*X2-i P2/IP2-P1) ) * { X2-X 1 ) 

P = M + 1 
GO TO 6C 

685 P=SIG1 C J 1+SIG21 J) 

IF( P)686,688,688 

686 I F ( 10 . 0 **{- 8 . 0 )-AeSF(P) )6e7, 693,693 

687 X1=X{J) 

P1=P 

X(J)=X2-(P2/(P 2-P 1 1 )■ (X2-X1 1 
GO TO 60 

688 I F C 10.0**1-8.01 -A BSF(P) 1689,693,693 

689 X2 = X ( J ) 

P2=P 

X( J ) = X2- ( P2/C P2-P 1 1 )*(X2-X1 ) 

GO TO 6C 
693 N-0 
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KTENS=1 



GO TO 698 

694 IF ( SI G 1 ( J ) +SIG2 ( J )) 695 , 698, 698 

695 IF( C. 5«BETA*TC-3. 141593)6955,696,696 
6955 KTENS=0 

GO TO 675 

696 KONTACT = C 
GO TO 675 

6 98000 V ( J )=C. 04 10 639# ( SLNESS**2 ) » ( ( S IG1 ( J ) *SIG2 ( J ) )* ( VC* 1 . 0 ) - 
1 (SIG1 ( J )- S IG2 ( J ) ) /ECC) 

AK ( J , 1)=DELT*( VEL+CI J) •DELT»CCV( J) > 

70 AK( J,2 )=CELT*DDVI J) 

DISPLt I )=X( 1 ) 

FORCE! I ) = .0410639«( SLNESS»*2 )•( STRS1< I )*STRS2 ( I) ) 

VEL=VEL* ( AK{ 1 , 2 )♦ 2. 0 • AK { 2 , 2 ) *2. 0* AK ( 3 , 2 ) ♦ AK ( 4 , 2 ) )/6.0 
OEFL ( 1 + 1 )=DEFL l I) ♦ ( AK ( 1 , 1 ) *2. 0* AK ( 2 , 1 ) *2 .0*AK ( 3 , 1 ) *AK ( 4 , 1 >)/6.0 
IF ( 1-900)75,80,80 
75 TCI + 1 ) = T( I ) *DELT 
1=1*1 
GO TO 4 C 

80 WRITE OUTPUT TAPE 4,82 
82 FORMAT ( 1 H 1 ) 

WRITE OUTPUT TAPE 4 , 85 , SLNE S S , ALPHA , E CC , B ET A , A , S IG Y , C ELT 

» 

850F0RMAT ( 14H OUTPUT OATA //3X 1 2HSLENDERNESS-F 5 . 0 ,2 X6HA LPHA=F 5 . 2 , 

12X4HECC=F5.3,2X5H8ETA=F4.2,2X2HA=F5. 1 ,2X5HSIGY«F6. 3,2X5H0ELT =F5. 3 
// ) 



2 
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WRITE OUTPUT TAPE 4,90 



9CCFORMA T ( 4 X 4HT I ME, 4 X12HD IS PLACEMENT, 4X10H0E FLECTION, 4X8HSTRAIN 1 , 
16X8HSTRESS 1 , 6X8HSTRA IN 2.6X8HSTRESS 2, 6X5HFORCE// ) 

00 92 1=1 ,900 

92 DEFL ( I ) = ECC»DEFL( I ) 

DO 95 1*1, 900,15 

950WRITE OUTPUT TAPE 4 , 98, T ( I ) ,0 ISPL (I) , DEFL ( 1 1 , STRNl ( I ) , STRS 1 ( I ) , 
1 STRN2( I ) ,STRS2( I) , FORCE! I ) 

98 FORMAT (3XF6.2»7E1 4.4) 

CALL GRAPH1900.T, STRM.8 ) 

CALL GRAPH(900, T, STRN2,8) 

CALL GRAPH(9C0,T, STRS1,8 ) 

CALL GRAPH(900,T, STRS2.8) 

CALL GRAPH! 900, T, 0EFL,8) 

CALL GRAPH(900,T, FORCE, 8) 

STOP 

END 

SUBROLTINE ST RE SS ( E PS I , E PS I A , S I GA , S I G ) 

COMMON A , S I GY 
I F { EPSI 120,20,21 

20 SIG=SIGA4(EPSI-EPSIA) 

GO TO 33 

21 IFl EPSI-EPSIAJ22, 22,23 

22 S1G=SIGA4(EPSI-EPSIA) 
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GO TO 33 



23 SIGB=$IGA + (EPSI-*EPSIA) 

8=3.0/(7.0*(SIGY**(A-1.0) )) 

S IG 1 =EPS I 

« 

2U SIG2 = SIG1-(SIGUB*(SIG1**A)-EPSI ) / ( A*B* < S IG 1 * * ( A-l .0) ) ♦ 1 .0 ) 
IF(ABSF(SIG2*(10. 0**{-6.0) ) >-ABSF ( S IG2-S IG 1 ) >29,30,30 

29 SIG1 =S IG2 
GO TO 24 

30 IF (SIG2-SIG8)31,3 1,32 

31 S I G=S IG2 
GO TO 33 

32 SIG=SIGB 

33 RETURN 
ENO 
END 
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d. Notation for Program COL 2 

In the notation for Program COL 2 the subscript J re- 
fers to quantities within the Runge-Kutta prediction loop 
which are calculated for intermediate times within the time 
increment AT. The variable quantities ending in B corre- 
spond to a time T - AT and the variable quantities ending 
in A correspond to a time T. The notation for Program COL 2 
is as follows: 



Computer program no t ation f or P rogram COL 2 



Computer 

Program 

Term 


Corresponding 
terms in table 
of symbols 


Computer 

Program 

Term 


Corresponding 
terms in table 
of symbols 


A 


n 


SIG1( J) 


s, 


ALPHA 


cx 


SIG2( J) 


Sa 


BETA 


$ 


SLNESS 


L/^ 


DDV ( J) 


i « 

V 


STRN1A or B 




DEPLA or B 


V 


STRN2A or 3 




DEFLP 


yr/<^ 


ST RSI A or B 


s, 


DELT 


AT 


STRS2A or 3 


s z 


DISPLA or B 


x/x £ 


STRS1M 


Max S, (tensile 
or compressive) 


ECC 




STRS2M 


Max S z (tensile 
or compressive 


EPSI1( J) 


6, 


T 


T 


EPSI2( J) 




VELA or B 


V 


PM AX 




X(J) 


x/x E 


SI GY 


Sy 
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e. Operating instructions for Program COL 2 

To operate Program COL 2 use the cards Including Sub- 
routine STRESS in the order shown in the listing followed by 
one blank card, then the input data card. 

The read-write statements in the listing are written for 
magnetic tape Input and output as follows: 

(1) Input (READ INPUT TAPE 2) 

The input to the program is contained on one card (7E10.3) 
and includes the following seven quantities in order: 

Input quant ity Format 

SLNESS(LA> ) E10.3 

ALPHA ( ) E10.3 

ECC (e/^) E10.3 



BETA ( Q ) 
A (n) 



E10.3 



E10.3 



SIGY (Sy) 
DELT (AT) 



£10.3 



E10.3 



(2) Output (WRITE OUTPUT TAPE 4) 



The output includes the input quatities and one value 
of each of the variable quantities including the time 



T of occurrence of the value as follows: 



Output quantity 
SLNE3S (L/^) 



Format 



F5.0 



ALPHA ( cx. ) 
ECC (e/^ ) 
BETA ( § ) 

A (n) 



F5.2 



E14.4 



F4.2 



F5.1 
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Output quantity 


Format 


SIGY (Sy) 


F6.3 


DELT (AT) 


F5.3 


PMAX (P m /P t ) 
T1 


E14.4 

F6.2 


STRS1M (S lfna J 
T2 


E14.4 

F6.2 


STRS2M (S lm J 
T3 


E14.4 

F6.2 


DEFLP (yr/o) 
T4 ' 


E14.4 

F6.2 
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PROGRAM COL 2 



C PROGRAM TO DETERMINE RESIDUAL DEFLECTION AND MAXIMUM FORCE 
C IN INELASTIC COLUMN RESULTING FROM HALF-SINE RAM DISPLACEMENT 
C PULSE APPLIED AT END WITH SMALL ECCENTRICITY TC INITIALLY 
C STRAIGHT COLUMN. WHEN RAM CAUSES TENSION COLUMN EXHIBITS FREE 
C VIBRATION. 

ODIMENSICN XU) „OD V ( 4 ) , EPS II C 4 ) , EPS 121 4 ) , S IG 1 ( 4 ) , SI G2 ( 4 > , AK ( 4 ,2 ) , 
1CU> 

COMMON A , S I GY 

READ INPUT TAPE 2 , 5 , SLNESS, AL PHA , ECC, BET A , A , S I C-Y ,D ELT 
5 FORMATI7EI0.3) 

1C T=0.0 

Cl SPL A=0 . 0 
VEL A=0 .0 
DEFLA=C.C 
STRN1 e-o.o 
STRN2 E=0 . 0 
STRS1 e=c.c 
STRS2B=0.0 
PMAX=C . 0 
STRS1 M=C . C 
STRS2 M = C • 0 

i 

STRS2E=0 . 0 
WORK = C .0 
C(1 )=C.C 
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C ( 2 ) =C .5 



C(3)*0.5 
C ( 4 ) = 1.0 
KENERG® 1 
KT ENS= 1 
KONT ACT = 1 
M= 0 

40 DO 70 J=1 ,4 
TC=T+C ( J ) *DELT 
VC=DEFLA*Cl J)*AK( J-1, 1 ) 

IF(0.5*BETA*TC-3. 141593)50,55,55 

50 X( J)=ALPhA*SINF(0.5*8ETA*TC) 

GO TO 60 

55 X(J)=C.C 

600EPSI 1 (J)=12.1761/ ( SLNESS**2 > * ( X( J )-0. 5* ( ECC**2 )• VC* *2* 
1ECC*VC*( 1 .O-ECC)) 

I F (J— 1 )62,62,63 

62 CALL STRESS(EPSI1 ( J ) , STRN IB, STRS 1 B,SI G1 ( J ) ) 

STRN1 A = EPSI 1 ( 1 ) 

STRS1 A=SIG1< 1 ) 

GO TO 64 

63 CALL STRESStEPSIl ( J ) , STRN 1 A, STRS 1 A, SI GU J) ) 

640EPSI2 l J ) = 12. 1761/ ( SLNESS**2 )* ( X ( J )-0. 5* ( ECC**2 )* VC* *2- 
. 1 ECC *VC* ( 1 .O + ECC)) 

IF ( J- 1 ).65 , 65, 66 
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65 CALL STRESS(EPSI2 ( J) , STRN2B, STRS2B,SIG2< J) ) 



STRN2 A = E PS I 2 ( 1 ) 

STRS2 A = S I G2 l 1 ) 

GO TO 661 

66 CALL STRESS<EPSI2(J),STRN2A,STRS2A,SIG2(J>> 
661 IFIKTENS-l 167,665 ,665 

665 IF(KOKTACT— 1 167 ,694,694 

67 IF(M-1)675,68,685 
675 X1=X(J) 

P1=SIG1(J)+SIG2CJ ) 

X( J)=CISPLB 
M = M+1 
GO TO 6C 

68 X2=X(J) 

P 2 = S I C- 1 ( J1+SIG21J ) 

XIJ) =X2- ( P2/( P2-P 1 ) ) «(X2-X1 ) 

M=M+1 
GO TO 60 

685 P = SIG1 ( J) +SIG2U) 

I F ( P ) 6 86 * 688 , 688 

686 I F ( 1 0 . C*» 1-8.0) -A 8SFi P) } 6 87,693,695 

687 Xl=X(J) 

P1 = P 

XU) = X2-(P2/1P2-P i ) 1 * ( X2-X 1 J 
GO TO 60 
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688 I F ( 10.0«*(-*8.0)-A8SF(P) )689,693»693 

689 X2=X(J) 

P2 = P 

X( J ) = X2- ( P2/( P2-P 1 ) ) * ( X2-X1 ) 

GO TO 6C 

693 M=0 
KTENS= 1 
GO TO 698 

694 IF(SIG1(J)+SIG2(J ) >695,698,698 

695 IF (0. 5*BETA*TC-3. 141593)6955,696,696 
6955 KTENS=0 

GO TO 675 

696 KONT ACT = C 
GO TO 675 

6980DDV<J)=0. 04 10639* ( SLNESS**2 )•( (SIG1 (J)+S«IG2tJ> 1MVC+1.C) 

I 

1 (SIG1 (J)-SIG2( J)> /ECC> 

AK( J, 1 )=CELT* ( VEL A + C ( J)»DELT*0DV( J ) ) 

70 AK(J,2)=DELT*DDVl J) 

IF (0.5*BETA*T-3.1 41593)702,704,704 
7020^0RK=fc0RK+0. 5*( < S TRS 1 A+S TRS2 A ) ♦ ( STRS 1 B+S TRS2B ) )* 

1 (ALPHA*(SINF(0.5*BETA*T)-SINF(0.5*BETA*(T-OELT) ) )) 

704 I F ( PMAX-C. 0410639* ( SLNESS**2 ) * I STRS 1 A + STR S2 A ) ) 7 1 , 7 2 , 7 2 

71 PMAX^C .04 10639* ( SLNESS**2 ) * { STRS 1 A+STRS2A) 

T 1 = T 

72 I F ( SIRS) V-STRS 1A) 73,74,74 
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73 

74 

75 

751 

752 

755 

756 

757 

758 

759 

76 

77 

78 

79 

80 



STRSlf* = STRS1A 



STRS1 E=STRS1M 
T2 = T 

IF(ABSF(STRS2M)-ABSF{STRS2A) >75,751,751 
STRS2N = ST RS2A 
T3 = T 

I F( STRS2E-STRS2A) 752,755,755 
STRS2 E = ST RS2A 
IF(KENERG-1 >756,76,76 
IF(VEIB)758» 759,7 57 
IF(VELA)759,7 59,8.1 
IF(VELA)81,759,759 
B=3.0/(7.0*( S IGY* • ( A- 1 .0 ) l ) 

ENERG 1= ( A*B/ (A + l .0) ) *STR$ 1 E** ( A+ 1 . 0) *0 . 5*S7R S 1 B ** 2 
ENERG2= ( A*8 / ( A +1 .0) ) *ST RS2E* * ( A + 1 .0 ) +0. 5*STRS2B**2 
ENERGT=0. 0821278* ( SLNESS**2 > *< ENERG 1 +ENERG2 ) 

GO TO 84 

IF ( KONTACT- 1 >77 ,8 1 ,81 
BENDM1=STRS1B-STRS28 
BEN0M2 = STRS1 A-STR S2 A 
I F ( STRS1B >78,78,79 
IF(STRS1A)81,80,80 
IF ( STRS 1A >80, 80,8 1 

DEFLP= ( DEFL A -BEN CN2 / ( 8ENDM2- BENDM 1 ) * ( OEFL A-OEF L B ) )*ECC 
T4 = T 
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KENERG=0 



81 T=T+DELT 

0 I SPLE=D I SPLA 
DISPLA=X(1) 

VELB=VEL A 

VEL A= VELB+ ( AK ( 1 ,2)+2.0«AK(2,2)+2.0*AK(3,2)+AK(4,2) >/6.0 
DEFLB=0EFLA 

DEFLA = DEFLB ♦( AMI , 1 ) +2 .0*AKC 2, 1 )*2. 0*AK ( 3, 1 ) +AK { 4, 1 ) ) / 6 . C 

STRN 1 B = STRN 1 A 

STRN2 E s STRN2A 

STRS 1 E=STRS1 A 

STRS2B=STRS2A 

GO TO 40 

84 WRITE OUTPUT TAPE 4,85,SLNESS, tCC , BETA ,A , $1 GY ,0E LT 
850F0RMAT ( / // 1 4H OUTPUT OATA //12HSLEN0ERNESS«F5.0,2X6H ECC=F5.2, 
12X5HBETA=F4.2,2X2HA=F5. 1 ,2X5H$IGY*F6. 3, 2X5HOELT*F5 . 3/ / ) 

WRITE OUTPUT TAPE 4,86 

860 F0RMATI6X5H ALPHA, 7X1 1 HOE FORMAT ION, 3X9HMAX FORCF, 5X 1 2HMAX STRESS 1, 
12X12HNAX STRESS 2//) 

87 WRITE OUTPUT TAPE 4 ,88, ALPHA, OEFLP , PM AX, STRS1 M r $TRS2M 

88 F0RMATI5E 14.4) 

WRITE OUTPUT TAPE 4 , 885 , T4, T1 , T2, T3 . 

885 FORMAT (20X1HIF6. 2, 1 H ) , 6X1 H ( F6 .2 , 1H ) ,6X1 H ( F6 . 2 , 1H ) , 6X 1 H ( F6. 2, 1H ) ) 

* 

WRITE OUTPUT TAPE 4 , 887, WCRK, ENERGT 
887 FORMAT { 25 X5HW0RK* E 1 2 .4, 1 0X7H ENERGY* E 1 2.4/) 
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STOP 



ENO 

SUBROUTINE STRESS (EPS I, EPSIA.S IGA, SIG) 

COMMON A, SIGY 
IF(EPSI)20»20»21 

20 SIG=S IGA+(EPSI-EPSIA> 

GO TO 33 

21 IF(EPSI-EPSIA)22, 22,23 

22 SIG=SIGA+(EPSI-EPSIA) 

GO TO 33 

23 SIGB=SIGA+(EPSI-EPSIA) 

B=3.0/(7.0*(SIGY*»(A-1 .0) ) ) 

S IG1 =EPS I 

24 SIG2 = SIGI-(SIGl+B*(SIG1**A)-EPS , n/(A«0«(Siqi*«(A-1 .01 + 1 . 0 ) 
IF(ABSF(SIG2*( 10. 0«* (-8.0)) )-ABSF( SIG2-SIG1 1)29,30,30 

29 SI G1=SIG2 
GO TO 24 

30 IF(SIG2-SIG8)31,31,32 

31 S I G=S IG2 
GO TO 33 

32 S IG*S IG8 

33 RETURN . . 

END 

END 
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APPENDIX V 



Maximum Force. Residual De flection and Maximum Stress Data 

Tabulated on the following pages are the maximum force, 
residual deflection, and maximum stress data for the three 
slenderness ratios and two eccentricity ratios for all values 
of CK and Q which were Investigated for pulse loading. The 
values of maximum lateral deflection that were available are 
also Included in the tabulation. These values were only avail- 
able from the solutions to computer Program COL 1 used to ob- 
tain graphical results. Thus, only a limited number of these 
values are Included in the tabulation. 

The dimensionless time Increments AT used for the numer- 
ical integration were the same for a given value of in 
all the solutions for both computer programs. These values 
are as follows: 



$ 


AT 


0.01 


0.10 


0.05 


0.10 


0.10 


0.08 


0.20 


0.04 


0.50 


0.02 


1.00 


0.01 



The values of AT are more than adequate to maintain the 
accuracy of the numerical scheme. The minimum number of 
increments for a natural period of first mode lateral vibra- 
tion of the column in the numerical process is 62. The values of 
AT were selected to display the complete history of the 
dynamic response to the pulse loading for the graphs shown 
in App. Ill using the maximum 900 points permitted in Sub- 
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routine GRAPH * 

The values of maximum stress tabulated on the following 
pages are maximum compressive stresses if there Is no sign 
before the number, or maximum tensile stresses if there Is 
a negative sigh before the number. In either case they are 
the maximum absolute values of stress that occurred In Spring 
1 or Spring 2 of the column model. 
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Maximum Force, Residual Deflection and Maximum Stress Data 



2 I/^- 1Q0 e/^« 0,05 



$ 


Maximum 

Force 

(p„/p 6 ) 


Residual 

Deflection 

(y*-/?) 


Maximum 

Deflection 

(yn/^) 


Stress 1 
Maximum 


Stress 2 
Maximum 
< 


0.01 


0.953 


2.12x10"* 




0.00282 


0.00074 


0.05 


1.001 


2.24x10”* 


— • 


0.00283 


0.00078 


0.10 


1.121 


2.50x10'* 


— 


0.00287 


0.00100 


0.20 


1.319 


2.77x10-* 


— 


0.00289 


- 0.00129 


0.50 


1.65$ 


2.40x10"* 




0.00286 


- 0.00233 


1.00 


1.935 


1.69x10"* 


— - ' 


0.0027$ 


0.00213 



<=<- 2.5 




e/^a 0.05 



$ 


Maximum 
Force 
<*U/fV ) 


Residual 

Deflection 


Maximum 

Deflection 

(y»n/£) 


Stress 1 
Maximum 


Stress 2 
■Maximum 
( s 2.max) 


0.01 


0.953 


5.89x10"* 


~ 


0.00312 


0.00084 


0.05 


1.029 


6.45x10"* 


— - 


0.00315 


-0.00099 


0.10 


1.189 


7.66x10"* 


— 


0.00320 


-0.00132 


0.20 


1.470 


9.41x10"* 


— 


0.00327 


-0.00188 


0.50 


1.958 


8.88x10”* 


_ — i 


0.00325 


r o . 00275 


1.00 


2.331 


9.00x10"*’ 




0.00329 


-0.00263 
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Maximum Force, Residual Deflection and Maximum Stress Data 



«<- 3 . ' L/^>- 100 e/^» 0.05 



$ 


Maximum 
Force 
(Pm/P c ) 


Residual 
Deflection 
(yr/^) . 


Maximum 

Deflection 

(ym/9) 


Stress 1 
Maximum 

(^irnax) 


Stress 2 
Maximum 


0.01 


0.954 


1.14x10"’ 


--- 


0.00334 


-0.00110 


0.05 


1.069 


1.28x10"' 


- — — 


0.00337 


-0.00130 


0.10 


1.252 


1.54x10-' 


— — — 


0.00344 


-0.00167 


0.20 


1.602 


1.98x10“ 


— — _ 


0.00352 


-0.00236 


0.50 


2.237- 


2.20x10“ 


— — » 


0.00356 


-0.00310 


1.00 


2.633 


2.05x10“ 


— : 


0.00360 


-0.00313 



«*>- 3.5 




e/^a 0.05 



> 


Maximum 

Force 

(p„/r, ) 


Residual 

Deflection 

(yr/^) 


Maximum 

Deflection 

(y»n/^) 


Stress 1 
Maximum 


Stress 2 
-Maximum 

( ax) 


0.01 


0.955 


1.81x10“ 




0.00349 


-0.00131 


0.05 


1.110 


2.05x10“ 




0.00354 


-0.00155 


0.10 


1.314 


2.47x10“ 


_ — _ 


0.00360 


-0.00195 


0.20 


1.724 


3.25x10“ 


— — - 


0.00370 


-0.00268 


0.50 


2.447 


3.83x10“ 




0.00377 


-0.00340 


1.00 


2.851 


2.57x10“ 




0.00380 


-0.00352 
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Maximum Force, Residual Deflection and Maximum Stress Data 



c <- 4 I/^-100 e/^„ 0.05 



$ 


Maximum 

Force 

(pVp e ) 


Residual 

Deflection 

(yr/9) 


Maximum 

Deflection 

(ym/9) 


Stress 1 
Maximum 

( S Unnax) 


Stress 2 
Maximum 

( nax) 


0.01 


0.956 


2.54x10"' 


■p— .. „ 


0.00361 


-0.00149 


0.05 


1.147 


2.39x10“' 


— 


0.00366 


-0.00175 


0.10 


1.362 


3.49x10“' 




0.00373 


-0.00219 


0.20 


1.336 


4.56x10"' 





0.00333 


-0.00231 


0.50 


2.617 


5.57x10“' 


— 


0.00392 


-0.00367 


1.00 


3.011 


2.93x10“' 




0.00394 


-0.00331 



o<» 1 70 e/^a 0.05 





Maximum 

Force 

(p„/p, ) 


Residual 

Deflection 

(yr/^) 


Maximum 

Deflection 

(y»n/^) 


Stress 1 
Maximum 
( 


Stress 2 
-Maximum 
( 6 2, max) 


0.01 


0.341 


3. 17x10" 1 


• 


0.00315 


0.00149 


0.05 


0.340 


3 . 19xl0 -i 


— 


0.00315 


0.00150 


0.10 


0.347 


3. 13x10“ z 


— 


0.00315 


0.00153 


0.20 


0.363 


4. 13x10“ 1 


— 


0.00324 


0.00163 


0.50 


0.934 


.2.33x10"* 




0.00312 


0.00190 


1.00 


0.932 


6.24x10"’ 


— 


0.00272 


0.00224 



V-5 



APPENDIX V 



Maximum Force, Residual Deflection and Maximum Stress Data 



1.5 



70 



e/^. 0.05 



$ 


Maximum 
Force 
(Pm/P 6 ) 


Residual 

Deflection 

(y*-/Q) • 


Maximum 

Deflection 


Stress 1 
Maximum 


Stress 2 
Maximum 

( s zma^ 


0.01 


0.841 


2.10x10“' 





0.00381 


0.00149 


0.05 


0.854 


2.09x10"' 


— 


0.00381 


0.00153 


0.10 


0.898 


2.10x10“' 




0.00381 


0.00164 


0.20 


0.987 


2.10x10“' 





0.00381 


0.00184 


0.50 


1.206 


2.04x10“' 


— 


0.00381 


-0.00260 


1.00 


1.337 


1 . 03 x 10 “' 


- — ■ 


0.00364 


0.00314 



2 Iy^ 70 e/^= 0.05 





Maximum 

Force 

) 


Residual 
Deflection 
( yr/^>) 


Maximum 

Deflection 

(y*n/0 


Stress 1 
Maximum 


Stress 2 
•Maximum 
( s 2.max) 


0.01 


0.841 


3.79x10“' 


- — 


0.00404 


0.00149 


0.05 


0.878 


3.86x10“' 


* * 


0.00405 


0.00.158 


0.10 


0.954 


3.96x10“' 


— 


0.00406 


0.00176 


0.20 


1.080 


4.05x10“' 


. 


0.00407 


-0.00228 


0.50 


1.377 


3.43x10“' 




0.00403 


-0.00368 


1.00 


1.516 


2.22x10“' 


— 


0.00402 


0.00364 
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Maximum Force, Residual Deflection and Maximum Stress Data 



c<- 2.5 L/^>- 70 e/^„ 0.05 



§ 


Maximum 
Force 1 

(p~,/p t ) 


Residual . 
Deflection 
(y*-/?) 


Maximum 

Deflection 

(ym/^) 


Stress 1 
Maximum 


Stress 2 
Maximum 


0.01 


0.641 


5.36x10“' 


— 


0.00416 


0.00149 


0.05 


0.697 


3.52x10“' . 


— 


0.00419 


0.00161 


0.10 


0.995 


5.77x10“' 





0.00421 


-0.00202 


0.20 


1.167 


. 5.99x10“' 


1.67 


0.00423 


-o„ 00306 


0.50 


1.469 


5.03x10*' 


— 


0.00422 


-0.00413 


1.00 


1.619 


3.06x10“' 


• 


0.00422 


0.00393 



1 , ly^a 50 6/^a 0.05 





Maximum 

Force 

UU/r. ) 


Residual 

Deflection 

(yr/^). 


Maximum 

Deflection 

(y»n/0 


Stress 1 
Maximum 
( 


Stress 2 
■Maximum 

( SlLmax) 


0.01 


0.565 


2.67x10“' 


— 


0.00417 


0.00234 


0.05 


0.592 


2.69x10“' 


— 


0.00416 


0.00240 


0.10 


0.614 


2.70x10 h 


— 


0.00419 


0.00251 


0.20 


0.655 


2.63x10“' 


— 


0.00417 


0.00275 


0.50 


0.731 • 


2.52x10“' 


— 


0.00419 


0.00329 


1.00 


0.774 


6.60x10“* 




0.00396 


0.00366 
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Maximum Force, Residual Deflection and Maximum Stress Data 
c<- 1.5 1/^m 50 e/^» 0.05 



$ 


Maximum 

Force 

(p,yp s ) 


Residual 

Deflection 

(y»-/Q) 


Maximum 

Deflection 

(ym/^) 


Stress 1 
Maximum 

( ^imax) 


Stress 2 
Maximum 


0.01 


0.535 


5.17x10*' 


— 


0.00446 


0.00234 


0.05 


0.612 


5 . 16x10“' 




0.00446 


0.00251 


0.10 


0.654 


5. 13x10“' 




0.00446 


0.00275 


0.20 


0.708 


5 . 05 x 10 ”' 


. 


0.00446 


0.00307 


0.50 


0.822 


4.03x10”' 


— 


0.00445 


0.00383 


1.00 


0.864 


1.60x10“' 


— - 


0.00435 


0.00413 



o<» 2 1^5- 50 e/^» 0.05 





Maximum 
Force 
(?„/!>« ) 


Residual 

Deflection 

(yr/^) 


Maximum 

Deflection 

(y*n/£) 


Stress 1 
Maximum 

( ®imax) 


Stress 2 
Maximum 
( ^ 2. max) 


0.01 


0.535 


7 . 27 x 10 "' 


— 


0.00461 


0.00235 


0.05 


0.626 


7.34x10"' 


— 


0.00462 


0.00257 


0.10 


0.679 


7.37x10”' 


— 


0.00462 


0.00287 


0.20 


0.756 


7.18x10”' 


— 


0.00462 


0.00336 


0.50 


0.873 


4.93x10”' 




• 0.00461 


-0.00454 


1.00 


0.913 


2.07x10”' 


— 


0.00456 


0.00439 
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Maximum Force, Residual Deflection and Maximum Stress Data 



c<- 2.5 I/^- 50 e/^« o.05 



$ 


Maximum 
Force 
(Pm/Pc ) 


Residual 

Deflection 


Maximum 

Deflection 

(y.n/9) 


Stress 1 
Maximum 

(^irnax) 


Stress 2 
Maximum 


0.01 


0.586 


9.14x10”' 


— 


0.00472 


0.00236 


0.05 


0.645 


9.30x10”' 


— 


0.00473 


0.00269 


0.10 


0.700 


9.40x10“' 


— 


0.00473 


0.00300 


0.20 


0.794 


8.98x10“' 


— 


0.00474 


-0.00413 


0.50 


0.909 


6.17x10”' 


— 


’ 0.00475 


-0.00475 


1.00 


0.948 


2.41x10"' 




0.00471 


0.00456 



c<- 2 Iy^- 100 e/^a 0.01 





Maximum 

Force 

(p„/p« ) 


Residual 

Deflection 

(yr/^) 


Maximum 

Deflection 

(y»n/£) 


Stress 1 
Maximum 
( ) 


Stress 2 
•Maximum 

( ^ 2. max) 


0.01 


1.003 


2. 5 6x10” 1 


— 


0.00237 


0-. 00100 


0.05 


1.192 


3.03x10"* 


— 


0.00292 


0.00113 


0.10 


1.376 


3. 2 9x10“ 1 


1.64 


0.00295 


0.00149 


0.20 


. 1.616 


3.30x10"* 


1.81 


0.00295 


-0.00180 


0.50 


1.911 


3.07x10"* 


1.93 


0.00294 


-0.00231 


1.00 


1.988 . 


• 2. 95x10" 3 


0.69 


0.00250 


0.00235 



/ 
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Maximum Force, Residual Deflection and Maximum Stress Data 



2.5 • I /^- 100 0.01 



$ 


Maximum 
Force 
(Pm/P s ) 


Residual 

Deflection 

(yr/9). 


Maximum 

Deflection 


Stress 1 
Maximum 


Stress 2 
Maximum ' 


0.01 


1.017 


6.93x10“* 


— 


0.00317 


0.00100 


0.05 


1.258 


8.89x10"* 


— 


0.00325 


- 0.00139 


0.10 


1.486 


1.04x10"' 


— 


0.00331 


-0.00185 


0.20 


1.813 


1.17x10"' 


--- 


0.00334 


- 0.00250 


0.50 


2.264 


1.11x10"' 


— — — 


0.00335 


-0.00286 


1.00 


2.423 


2.29x10"* 




0.00306 


0.00287 



c<« 3 




e/^a 0.01 





Maximum 

Force 

(p„/p 5 ) 


Residual 

Deflection 

(Yr/^) 


Maximum 

Deflection 

(y*n /0 


Stress 1 
Maximum 
( 


Stress 2 
■ Maximum 
( 2 2 . max) 


0.01 


I.O 32 


1 . 30 x 10 "' 





0.00338 


-0.00115 


0.05 


1.320 


1 . 72 x 10 "' 


’ 


0.00348 


-0.00173 


0.10 


1.578 


2.07x10"' 





0.00354 


-0.00224 


0.20 


1.971 


2.47x10"' 


— 


0.00360 


-0.00274 


0.50 


2.535 


2.11x10 h 




0.00359 


-0.00326 


1.00 


2.744 


7.27x10"*- 


-— 


0.00346 


0.00327 
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Maximum Force, 

* 

<=*- 3.5 



$ 


Maximum 
Force ' 

(pVp s ) 


Residual 

Deflection 

(yr/ 9 ) 


Maximum 

Deflection 

(ym/ 9 ) 


Stress 1 
Maximum 

( S ur*ax) 


Stress 2 
Maximum 

( ^irnax^ 


0.01 


1.04S 


2. 00x10" ‘ 


— 


0.00353 


-0.00135 


0.05 


i.3ao 


2.69x10"' 


— 


0.00363 


-0.00199 


0.10 


1.666 


3.26x10“' 


: 


0.00370 


-0.00255 


0.20 


2.106 


3 .a 6 xlO“' 


— 


0.00377 


- 0.00292 


0.50 


2.739 


3.00x10“' 


- — 


0.00376 


-0.00356 


1.00 


2.961 


1.35x10“' 


— . 


0. 00372 


0.00354 



o-C. • 4 100 e/^a 0 • 01 



?■ 


Maximum 

Force 

( V*. ) 


Residual 

Deflection 

(yv/^) 


Maximum 

Deflection 


Stress 1 
Maximum 


Stress 2 
- Maximum 

( S^max) 


0.01 


1.062 


2.74x10“' 


--- 


0.00364 


-0.00151 


0.05 


1.433 


3.70x10“' 


— 


0.00375 


-0.00220 


0.10 


1.751 


4 . 52 x 10 “' 


2.67 


0.00333 


-0.00273 


0.20 


2.230 


5.49x10“' 


2.97 


0.00390 


-0.00303 


0.50 


2. 394 


3.94x10“' 


3.54 


0.0033a 


-0.00333 


1.00 


3.114 


1.93x10“' 


2.70 


0.00390 


0.00374 



APPENDIX V 

Residual Deflection and Maximum Stress Data 
l/<£- 100 e/^» 0.01 



V-ll 



APPENDIX V 



Maximum Force, Residual Deflection and Maximum Stress Data 



• 1 70 0,01 



$ 


Maximum 
Force 
<P~/P t ) . 


Residual 

Deflection 

(yr/9) 


Maximum 

Deflection 

(y.n/9) 


Stress 1 
Maximum 


Stress 2 
Maximum 


0.01 


0.916 


2.54xl0‘* 


- — 


0.00309 


0.00197 


0.05 


0.927 


2. 89x10- z 


— 


0.00313 


0.00201 


0.10 


0.947 


3.63x10”* 


0.49 


0.00320 


0.00206 


0.20 


0.973 


2.28x10”* 


0.45 


0.00306 


0.00217 


0.50 


0.991 


3. 94x1 0” 3 


0.20 


0.00264 


0.00232 


1.00 


0.994 


Iil3xl0" 3 


0.08 


0.00252 


0.00242 



<=<- 1.5 



L/^ 70 



a/^» 0.01 





Maximum 

Force 

dUA. ) 


Residual 

Deflection 

(yr/^) 


Maximum 

Deflection 

(y*n/^) 


Stress 1 
Maximum 

(^irnax) 


Stress 2 
. Maximum 
(Sz.max) 


0.01 


0.922 


2.24x10”' 


— 


0.00383 


0.00199 


0.05 


0.994 


2.24x10”' 


— 


0.00383 


0.00212 


0.10 


1.072 


2.22x10"' 


— 


0.00383 


0.00230 


0.20 


1.181 


2.17x10“' 


— 


0.00383 


0.00261 


0.50 


1.333 


1.60x10”' 





0.00377 


0.00315 


1.00 


1.364 


2.11x10”* 


— 


0.00345 


0.00334 
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Maximum Force, Residual Deflection and Maximum Stress Data 



2 



1 /^- 70 



e/^Ja 0.01 



$ 


Maximum 

Force 

(P m /P s ) 


Residual 

Deflection 

(yr/^) 


Maximum 

Deflection 

(ym/9) 


Stress 1 
Maximum 

( s irnax) 


Stress 2 
Maximum 


0.01 


0.930 


4.01x10“' 


mm mm mm 


0.00406 


0.00200 


0.05 


1.041 


4 . 17 x 10 "' 


— 


0.00408 


0.00222 


0.10 


1.151 


4.17x10“' 


— 


0.00408 


0.00250 


0.20 


1.296 


4 . 02 x 10 “' 


— 


0.00408 


- 0.00309 


0.50 


1.497 


3.04x10“' 


— . 


0.00408 


0.00359 


1.00 


1.542 


5.20xl0“ i 


— 


0.00388 


0.00380 



c<- 2.5 Iv^« 7 Q e/^» 0.01 





Maximum 

Force 

(P m /r e ) 


Residual 

Deflection 

(y»-/^) 


Maximum 

Deflection 

(y*«/£) 


Stress 1 
Maximum 


Stress 2 
-Maximum 

( ^ z. max) 


0.01 


0.939 


5.62x10“' 




0.00420 


0.00201 


0.05 


1.079 


6.02x10“' 


• 


0.00423 


0.00231 


0.10 


1.210 


6.17x10“' 


1.85 


0.00424 


-0.00292 


0.20 


1.362 


6 . 03 x 10 “' 


1.95 


0.00425 


- 0.00364 


0.50 


1.594 


3 . 95 x 10 “' 


1.99 


0.00425 


-0.00396 


1.00 


1.642 


7 . 93 x 10 “ 2 - 


0.72 


0.00413 


0.00405 
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Maximum Force, Residual Deflection and Maximum Stress Data 



o<.- 1 



L/^ ■ 50 



e/^ a 0.01 



§ 


Maxi nun 
Force 

(?m/? S ) 


Residual 

Deflection 

(yr/9) 


Maximum 

Deflection 

(*»/?> 


Stress 1 
Maxi nun 

( ^lrnax^ 


Stress 2 
Maxi nun 

I 8 *™*) 


0.01 


0.633 


2.71x10"' 


— 


0.00419 


0.00234 


0.05 


0.658 


2.74x10"' 


— 


0.00420 


0.00297 


0.10 


0.689 


2.75x10“' 


0.67 


0.00421 


0.00313 


0.20 


0.730 


2.65x10“' 


0.75 


0.00421 


0.00337 


0.50 


0.775 


1.15x10“' 


0.48 


.0.00403 


0.00367 


1.00 


0.782 


1.78xlO“ l 


0.10 


0.00384 


0.00378 



o 1,5 Iy^» 50 e/^s 0,01 





Maxi nun 
Force 

(?„/?, ) 


Residual 

Deflection 

( >v/^>) 


Maximus 

Deflection 

(y»n/ 


Stress 1 
Maximum 

( ^irna.x^ 


Stress 2 
. Maximum 
( S^max) 


0.01 


0.636 


5.30x10"' 


— 


0.00447 


0.00286 


0.05 


0.637 


5.22x10“' 


— 


0.00447 


0.00311 


0.10 


0.735 


5.02x10“' 


1.12 


0.00447 


0.00338 


0.20 


0.791 


4.57x10“' 


1.10 


0.00447 


0.00369 


0.50 


0.859 


2.62x10"' 


0.90 


• 0.00443 


0.00411 


1.00 


0.870 


3 . 66x10“* 


0.16 


0.00426 


0.00422 
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Maximum Force, Residual Deflection and Maximum Stress Data 

■ 1 L/^ m 50 6/^ a 0*01 



$ 


Maximum 

Force 

ou/p« ) 


Residual 

Deflection 

(y*~/<?) 


' Maximum 
Deflection 

(Ynr,/^) 


Stress 1 
Maximum 


Stress 2 
Maximum 


0.01 


0.633 


2.71x10“' 


— 


0.00419 


0.00284 


0.05 


0.658 


2.74x10“' 




0.00420 


0.00297 


0.10 


0.689 


2.75x10“' 


0.67 


0.00421 


0.00313 


0.20 


0.730 


• 2.65x10“' 


0.75 • 


0.00421 


0.00337 


0.50 


0.775 


1.15x10-' 


0.48 


0.00403 


0.00367 


1.00 


0.782 


1* 78x10“ * 


0.10 


0.00384 


0.00378 



o<- 1.5 Iy^. 50 e/q= 0.01 





Maximum 

Force 

owa* ) 


Residual 

Deflection 

(yr/^) 


Maximum 

Deflection 

(y*n/£) 


Stress 1 
Maximum 


Stress 2 
-Maximum 
( s 2.max) 


0.01 


0.636 


5.30x10“' 


mm mm 


0.00447 


0.00286 


0.05 


0.687 


5. 2 2x10“ 1 


mm mm mm 


0.00447 


0.00311 


0.10 


0.735 


5.02x10-' 


1.12 


0.00447 


0.00338 


0.20 


0.791 


4.57x10“' 


1.10 


0.00447 


0.00369 


0.50 


0.859 


2.62x10"' 


0.90 


• 0.00443 


0.00411 


1.00 


0.870 


3 • 66x10“* 


0.16 


0.00426 


0.00422 
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Maximum Force, Residual Deflection and Maximum Stress Data 



OC- 2 L/^>- 50 °* 01 



$ 


Maximum 

Force 

(P m /P s ) 


Residual j 
Deflection 
(yr/^) 


Maximum 

Deflection 

(y.n/<$) 


Stress 1 
Maximum 
( S \ max ) 


Stress 2 
Maximum 


0.01 


0.641 


7.44x10"' 


— 


0.00463 


0.00269 


0.05 


0.709 


7.44x10"' 


— 


0.00463 


0.00323 


0.10 


0.766 


7.17x10"' 


— 


0.00464 


0.00354 


0.20 


0.632 


6 . 27 x 10 “' 


— 


0.00463 


0.00391 


0.50 


0.907 


3.44x10"' 




0.00463 


0.00435 


1.00 


0.919 


4 . 44x10 " \ 


— 


0.00449 


0.00446 



c*.- 2.5 50 ®/^ 0.01 





Maximum 

Force 

(p„/p« ) 


Residual 

Deflection 

(yr/^) 


Maximum 

Deflection 

(y*n /0 


Stress 1 
Maximum 
( 


Stress 2 
. Maximum 
( s tmax) 


0.01 


0.646 


9.35x10"' 


— 


0.00473 


0.00291 


0.05 


0.729 


9.45x10"' 


- — 


0.00475 


0.00334 


0.10 


0.769 . 


9.19x10"' 


1.65 


0.00476 


0.00367 


0.20 


0.663 


7.94x10"' 


1.60 


0.00476 


-0.00427 


0.50 


0.941 


4.03x10"' 


1.26 


0.00476 


0.00452 


1.00 


0.952 


5.25xlO“ L 


0.22 


0.00465 


0.00462 



V-15 



APPENDIX VI 



Evaluation of the Approximat ions and of th e Accuracy of 
the Numerical Scheme 

From a consideration of the maximum lateral deflections 
encountered in the solutions the approximations made in de- 
veloping the equations for the column model may be evaluated. 
The maximum value of© for all solutions was less than 0.08 
radian. Based on a maximum © = 0 o 08 radian the errors are 
as follows: 

(a) 0.01 per cent error in (L/2-e tan© ) as a result 
of neglecting e tan© 

(b) 0.2 per cent error in Er a ine=e(l -cos-©- 3 by con- 
sidering sin©=© and neglecting e( 1-cos-©) 

(c) 0.25 per cent error in tan© and 0.11 per cent 
error in sin© as a result of considering 
tan © = sin©=© 

(d) less than 0.06 per cent error in (1-cos©) as a 
result of considering l-cos©=£© 4 

From this the approximations which have been made in deriv- 
ing the equations for the column model are reasonable and 
produce insignificant error in the results. 

As an additional check for cumulative error a compar- 
ison of the ram work done on the column model with the net 
energy input to each of the springs (a combination of strain 
energy of plastic deformation and the existing elastic strain 
energy) is made in all computer solutions for the results 
shown in Figs. 10 thru 14. These agree within 1 per cent 
thus verifying the compatibility of the equations and c'neck- 
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ing the accuracy of the numerical scheme for the time Incre- 
ments chosen. Work-energy considerations are not Included 
In the computer program for obtaining graphical results but 
the dimensionless time increment AT used for the numerical 
integration in both programs for a given value of was the 
same. 
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